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Abstract 

Troposcatter  data  were  collected  with  a  wideband  measurement  system  for  a  period  of 
5  months  spanning  pans  of  the  winter,  spring  and  summer  seasons  of  1991.  The  system 
provided  multipath  delay  measurements  with  a  2.5  nsec  resolution.  Measurements  of 
received  signal  level,  Doppler  frequency  shift  and  multipath  delay  spread  were  obtained  on 
the  161  km  troposcatter  path  from  Prospect  Hill  in  Waltham,  Massachusetts  to  Mt.  Tug  in 
Lebanon,  New  Hampshire.  These  observations  are  the  latest  in  a  longer  measurement 
program  to  characterize  the  propagation  channel  for  tropospheric  scatter  communications 
systems  operating  at  frequencies  near  16  GHz  (Ku  band).  Simultaneous  observations  at 
5  GHz  (C  band)  were  made  during  the  last  month  of  the  observing  program  to  provide 
calibration  and  a  statistical  comparison  with  earlier  observations  made  using  a  different 
Ku  band  transmission  system. 

The  results  obtained  with  the  wideband  system  were  similar  to  those  collected  earlier. 
The  intent  of  the  measurement  program  with  the  wideband  system  was  to  make  multipath 
delay  spread  measurements  with  better  delay  resolution  than  had  been  obtained  before. 
Model  calculations  had  suggested  that  the  delay  spread  should  be  significantly  narrower 
than  previously  reported  and  narrower  than  assumed  in  the  design  of  troposcatter 
modulation  systems.  The  prior  observations  made  during  this  program  employed  a  channel 
measurement  system  with  an  80  nsec  resolution.  The  observations  with  that  system  were 
consistent  with  the  hypothesis  that  the  atmospheric  channel  had  a  multipath  delay  structure 
that  was  significantly  narrower  than  had  been  assumed  and  too  narrow  to  observe  with  the 
existing  equipment.  The  newer  system  had  sufficient  delay  resolution  to  allow  a 
comparison  between  the  observations  and  predictions.  The  results  show  that  the 
atmospheric  channel  has  a  wider  bandwidth  than  previously  assumed  (a  narrower  multipath 
delay). 

Statistically,  the  median  delay  spread  values  for  the  1991  spring  time  period,  which 
had  with  the  longest  run  of  data,  ranged  from  98  nsec  for  the  afternoon  hours  to  1 10  nsec 
during  the  evening  for  times  with  scattering  by  clear  air  turbulence.  These  values  are 
narrower  than  the  173  to  184  nsec  median  values  obtained  with  the  older  system  during  the 
spring  of  1990  for  the  same  time  periods  respectively.  They  are  also  narrower  than  the  150 
to  160  nsec  values  that  have  been  used  in  system  design. 
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1.  Introduction 

1 . 1  Program  objectives 

The  objectives  of  this  study  were  to  use  the  wideband  troposcatter  system  to 
characterize  the  channel  for  use  in  a  digital  communications  system  at  operating  at  a  earner 
frequency  near  16  GHz. 

1 . 2  Summary  of  Results 

The  analysis  of  the  observations  made  with  the  wideband  troposcatter  measurement 
system  is  presented  in  the  Appendix  to  this  final  report. 

Data  collection  for  this  phase  of  the  experiment  started  in  January  1991  and  continued 
through  the  June  1991  with  occasional  outages  due  interference  in  several  of  the  TV 
channels  caused  (apparendy)  by  the  first  IF  of  the  Ku  band  receiver.  In  total,  800  hours  of 
observations  were  collected  during  the  winter  period,  1308  in  the  spring  and  359  in  the 
summer.  C-band  measurements  were  started  in  the  latter  pan  of  the  spring  and  continued 
during  the  summer  period.  A  total  of  642  hours  of  simultaneous  C  and  Ku  band 
measurements  were  obtained. 


Analysis  of  a  Wideband  Troposcatter  Propagation  Channel  at  the  Ku  Band 


Magnus  Wennemyr 
Master  of  Science  Thesis 
October  1991 

Thayer  School  of  Engineering 
Dartmouth  College,  Hanover,  New  Hampshire 


ABSTRACT 

Tropospheric  scatter  is  a  conglomerate  term  describu  g  the  course -of-propaganon  alteration  of 
electromagnetic  wave  energy  due  to  weather  mechanisms,  allowing  for  transhorizonal  communication  in  the 
microwave  ranges.  This  study  tests  the  viability  of  using  the  higher  frequency  ranges  (the  test  frequency 
was  15.731  GHz)  for  this  type  of  communication.  At  the  present  time,  troposcatter  communication 
systems  are  limited  to  the  lower  frequencies  (up  to  5GHz),  due  to  concerns  of  attenuation  at  the  higher 
frequencies  during  rainy  conditions. 

Data  were  collected  on  a  wideband  system  (400  MB  it/sec)  over  a  period  of  five  months.  This  system 
provides  an  unprecedented  resolution  of  turbulence  in  the  lower  troposphere.  During  the  last  month,  data 
were  also  collected  simultaneously  with  a  reference  C  band  system  (4.95  GHz).  Characteristics  of  the 
signal  level,  Doppler  shift  and  multipath  delay  were  recorded,  and  hourly  median  values  of  these  quantities 
were  analyzed.  Diurnal  and  seasonal  changes  were  observed,  and  the  data  categorized  according  to  the 
dominant  scattering  mode  for  each  hour.  Three  modes  were  considered:  turbulent  volume  scatter,  rain 
scatter,  and  ducting. 

The  observations  made  on  the  161km  long  path  between  Prospect  Hill,  MA,  and  MtTug,  NH.show  that 
rain  is  beneficial  for  wave  propagation  at  this  high  frequency.  The  increased  forward  scatter  associated  with 
the  rain  offsets  the  attenuation  caused  by  scattering  away  from  the  path  between  the  transmitter  and  receiver. 
Turbulent,  thin-layer,  volume  scattering  dominated  the  propagation  process  90  percent  of  the  time.  As  a 
theory  of  propagation,  this  clear  air  phenomenon  was  consistent  with  the  results.  Diurnal  variations  were 
strongest  in  the  summer,  significantly  increasing  both  the  Ku  and  C  band  signal  levels  during  afternoon 
hours.  Seasonal  variations  were  evident  as  well.  In  particular,  an  important  observation  of  seasonal  change 
was  a  delay  spread  increase  of  50  percent  from  winter  to  summer. 
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Chapter  1 
Introduction 

1.1  Introduction 

The  purpose  of  this  study  was  to  test  the  viability  of  using  a  higher  frequency  range  for  Tropospheric 
Scatter  itroposcaiter)  Communication.  Although  ihe  frequency  range  in  quesuon  -  near  and  around  ihe  ,est 
frequency  of  15.731  GHz  •  is  already  in  use  for  other  types  of  commumcaaon  and  remote  sensing  devices,  it 
has  been  assumed  that  troposcatter  would  be  too  severely  attenuated  at  such  small  wavelengths,  especial !v 
in  rain  conditions,  for  reliable  use  in  a  radio  communications  service.  iThe  15.4  to  17.3  GH/.  ranee  is 
allocated  for  RADAR  services:  the  closest  fixed  service  allocation  is  14.8  to  15  25  GH7  )  If  ccrirusunc 
theories  were  to  prove  accurate,  however,  the  funders  of  the  project  would  be  able  to  use  the  new 
frequencies.  These  theories  which  rely  on  a  theoreucal  framework,  rather  than  empirical  one.  predict  that 
propagauon  of  frequencies  in  the  15  GHz  range  will  be  enhanced  during  rain;  this  would  be  the  result  of  an 
increase  in  scattering  during  rain,  counteracting  the  accompanying  anenoauon.  Antenna  gain  can  be 
maintained  ai  higher  frequencies  while  reducing  the  antenna  size  -  such  a  rcducuon  would  make  troposcatter 
systems  pouable.  Furthermore,  used  as  a  tool  for  ohservauon  of  meteorological  processes,  this  bisiatic 
RADAR  system  has  resolved  the  vertical  Doppler  and  scaoercr  strength  components  of  the  lower 
troposphere  to  an  unprecedented  degree. 

1.2  Transhorizonal  Microwave  Propagation  in  the  Troposphere 

The  troposphere  is  the  innermost  layer  of  the  atmosphere,  ranging  from  ground  level  up  to  10- 16 
kilometers  height  depending  on  the  geographic  location.  The  mean  temperature  decreases  with  height 
within  this  Layer.  Beyond  the  tropopause  ■  a  constant  temperature  layer  between  the  troposphere  and  the 
stratosphere  -  the  temperature  again  increases.  All  weather  conditions  manifest  themselves  in  the 
troposphere,  significantly  affecting  the  propagauon  of  those  radio  waves  with  wavelengths  on  ihe  order  of 
the  vanous  inhomogeneities  in  the  refracuve  index.  The  unpredictable  and  complex  nature  of  weather 
obscures  interpretation  of  these  effects.  Most  significant  of  the  regions  where  these  disconunumcs  o \ i ' t . 
are  strong  layers  of  clear-air  turbulence  and  volumes  containing  the  different  forms  of  precipitation 

Tropospheric  scattering  is  a  conglomerate  term  describing  the  course  of  propagauon  alteration  of 
elcctromagneuc  wave  energy  due  to  these  weather  mechanisms,  allowing  for  transhorizonal  communicuticnv 
in  the  microwave  ranges.  By  several  mechanisms,  weather  •  defined  as  the  temperature,  pressure,  humidity 
wind  velocity,  and  precipitation  configuraaon  of  the  troposphere  -  enables  the  establishment  of 
communication  links  of  over  1000  km.  pathlengths.  Conversely,  by  the  same  phenomena  unpredictable 
changes  in  the  weather  can  also  inadvertantly  cause  interference  between  distant  co-channel  sen.  ices. 


Although  the  mtroducuon  of  satellite  communication  reduced  the  need  for  transhonzonal  tropo  scatter 
systems  tor  general  communications,  and  aside  from  their  use  in  special  circumstances  idue  to  the  low 
investment  costs  and  the  particular  capabilities  desirable  to  the  military),  they  are  useful  in  studying  the 
interference  p-oblem  for  systems  operating  at  the  same  range  of  frequencies  (e  g.  tne  European  Cammunuv 
COST  210  project  uses  tioposcatter  data  to  set  satellite  communtcauon  design  guidelines  for  avoiding 
interference  )  Moreover,  these  systems  are  used  for  the  study  of  meteorology;  by  idenufymg  the  processes 
that  alter  the  transmitted  signals  by  inference  from  differences  in  the  transmitted  and  received  signals, 
knowledge  is  gained  about  weather  processes  and  structure.  This  last  application  is  most  useful  in  dear  ac. 
during  such  conditions  the  different  components  of  the  multi-path  delay  signature  (i  c.  the  power  received 
a  function  of  the  transmiuo-receive  umelag)  can  rcahsucaJJy  be  considered  to  represent  a  verucal  profile  u 
the  seatierer  strength. 

1.3  The  Troposcatter  Link 


The  troposcaaer  link  shown  in  Fig  11  consisted  of  1)  a  transmitter  at  Prospect  Hill,  Waltham  MA, 
maintained  by  Rome  Laboratory  personnel  and  2)  a  receiving  stanon  on  MLTug.  Lebanon  NTt,  constructed 
and  maintained  by  Thayer  School  graduate  students  over  a  period  of  five  yean 


Figure  l.l.  Prospect  Hill  -  MLTug  path  profile  for  a  J/3  dfecuve  earth  radius.  Tlie  luwi-'.i 
beam  paths  shown  indicate  the  center  ut  the  conical  1/2  powei  Seuioriv  ,.-l  :k: 
uansmiiied  beams. 


The  angle  of  the  pointing  of  Lhe  antennas  with  the  horizon  is  small,  producing  a  very  small  common 
scattering  volume.  As  a  result,  this  system  is  ideal  both  for  studying  the  interference  phenomenon  which 
occurs  over  similar  paths  between  ground-based,  microwave  communication  systems,  and  for  achieving  the 
aforementioned  high  vertical  resoluuon.  The  armed  forces  use  this  type  of  system  because  cf  the  difficulty 
of  linking  with  the  small  scattering  volume  without  access  to  the  transmitter-receiver  sites,  providing  added 
secunry  from  interference  and  unwanted  listeners.  The  possibility  of  using  smaller  antennas  means  that 
these  systems  could  become  portable  by  vehicles. 

1.4  Project  Summary 

Building  of  the  link  was  begun  by  setting  up  the  microwave  (Ku  band  -  at  15.731  GHz)  receiver  [Levins 

1986] .  A  lower  frequency  receiver  (C  band  -  at  4  95  GHz)  was  set  up  as  a  standard  for  comparison  [Sull.vm 

1987] ,  This  frequency  was  chosen  to  be  high  enough  to  lie  within  the  same  subrange  of  turbulent  scales 
(the  inertial'  range  of  turbulence,  discussed  in  Section  2.3.2,  where  scatienng  strength  is  predicted  to  follow 
a  log-linear  frequency  dependence)  in  order  to  accurately  evaluate  clear  air  frequency  relationships.  The 
frequency  was  also  low  enough  so  that  a  substantial  amount  of  data  and  expenence  were  available  for 
comparison.  Since  the  C-band  had  only  this  purpose,  it  employed  simple,  continuous  wave  iCW) 
transmission.  To  provide  information  as  to  the  delay  characteristics  of  the  channel,  a  phase  modulated  1023 
bit  long  pseudo-random  .x>ise  (PRN)  code  at  12.5  MBit/s  (narrowband)  was  applied  to  the  experimental  Ku 
band.  [Suunboulis  1987  and  Anderson  1989  further  add  to  the  system  description.)  Later  a  400  MBit/s 
(wideband)  PRN  code  was  added  to  the  system  to  increase  the  spatial  resoluuor  (Chandra  1990). 

Preliminary  assessment  of  the  link  and  summary  of  the  narrowband  results  until  spring  1990  have  already 
been  made  [Hoppe  1989.  Prabhu  1991,  Crane  1991,  Crane&Prabhu  1991).  This  thesis  presents  an  analysis 
of  the  wideband  system  results. 

1.5  Personal  Contribution  to  Project 

A.side  from  the  daily  routine  of  collecting  and  analyzing  the  data,  only  minor  work  was  required  for 
maintenance  of  the  receiver  hardware.  The  wideband  receiver  installation  was  completed  dunng  the  fall  oi 
1990,  during  which  ume  the  transmitter  was  inoperauonal.  New  software  for  data  collecuon  was  developed, 
based  on  the  collection  program  for  the  narrow-band  system,  and  limitations  on  the  computer  at  the  receiver 
end  was  listed.  A  recurring  problem,  caused  by  a  faulty  clocking  of  Lhe  PRN  code  generator,  was 
eventually  fixed.  When  the  C  band  was  put  into  use  in  May.  interference  problems  interrupted  the 
collecuon  campaign.  Interference  from  the  C  band  onto  the  Ku  band  was  successfully  dealt  with,  by 
switching  the  systems  on  and  off  during  their  near-simultaneous  measurements  toward  the  end  of  die 
collecuon  campaign.  The  local  cable  TV  company  complained  of  interference,  for  which  reason  operation 
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was  shut  down  several  umw  while  a  solution  was  found.  The  data  processing  software  was  rewritten  m  ■/.> 
entirety,  and  a  method  greatly  diminishing  the  effects  of  the  unwanted,  irregular,  relative  drift  of  ihc 
frequency  Cesium  beam  standards  was  developed.  This  method  was  required  in  order  to  make  possible  a 
precise  defining  of  the  spatial  location  of  the  scatturers. 

1.6  Objective  of  the  Study 

The  purpose  of  this  report  is  to  present  the  effects  of  the  different  weather  phenomena  •  w  uh  an  emphasis 
on  ram  -  on  16  GHz  microwave  propagation  m  the  troposphere.  The  data  discussed  represents 
measurements  collected  during  a  six  month  period.  The  use  of  the  C  band  at  the  end  of  the  spring  was 
primarily  as  a  reference  signal.  Since  the  transmission  characteristics  depend  on  a  variety  of  factors,  one 
expects  deviation  from  theoretical  predictions,  wherefore  the  data  are  presented  in  a  stausncal  format,  and 
their  validity  investigated.  The  data  are  presented  in  a  format  to  be  used  for  a  communicauon  link.  High 
level  fields  are  in  that  case  of  little  interest;  rather,  low-level  fields  are  of  primary  concern.  The  most 
important  parameters  investigated  are  the  transmission  loss,  the  multipath  delay  spread  (placing  limitations 
on  channel  bandwidth),  and  the  Doppler  frequency  spread  (a  necessary  parameter  for  many  modulation 
schemes). 

1.7  Summary  of  Signal  Processing  and  Results 

In  order  to  characters.  "ommumcabon  channel  and  provide  the  means  for  interpretation  of  the 
meteorological  pheonomc^,  information  was  collected  of  the  power  received,  the  mean  Doppler  frequency 
shifts  and  the  Doppler  frequency  spread,  the  multipath  delay  spread,  as  well  as  power  and  Doppler 
measurements  at  different  delays.  Because  of  the  limitations  of  the  IBM  PC  used  for  data  collection  at  the 
reaver  site,  only  sixteen  scans  of  the  PRN  code  •  plus  data  for  a  spectrum  at  the  delay  with  the  maximum 
power  •  were  possible  during  an  hour.  The  in-phase  and  quadrature  components  of  the  power  were  sampled 
60  times  at  each  tap  of  the  code,  in  order  to  average  out  Rayleigh  fading  on  the  channel.  Hourly  medians  i 
the  different  quantities  were  then  arranged  in  various  ways  to  allow  investigation  of  the  diurnal  and  seasonal 
changes  in  the  data,  as  well  as  the  effect  due  to  different  weather  conditions  -  the  modes  of  propagation  In 
order  to  statistically  predict  the  occurrence  of  low-level  fields,  large  Doppler  spreads,  and  of  large  delay 
spreads,  distnbuuon  models  for  these  quanuues  were  fiued  to  the  data. 

1.7.1  Diurnal  and  Seasonal  Variations 

As  expected,  the  transmission  losses  were  smallest  during  the  afternoon  when  warm,  moist,  turbulent,  air 
i  "ibuted  more  to  the  forward  scattering  process.  The  vanauon  during  die  winter  was  minimal,  but 


increased  to  an  average  of  about  8  dB  during  the  summer  (June  1st- 20th).  There  was  a  progressive  increase 
in  this  vanauon  during  the  spring  months.  There  is  a  less  clearly  defined  diurnal  variauon  of  the  muitipath 
delay  spread;  it  tends  to  increase  dunng  the  later  hours  of  the  day,  following  die  increase  in  the  signal  level 
The  general  signal  level  increases  dunng  the  wanner  months  (a  difference  of  4-5  dB).  but  is  offset  by  a  very 
large  increase.  30-50%,  in  the  delay  spread  values.  The  algorithm  used  for  the  delay  spread  estimate  is  more 
sensitive  than  that  used  in  previous  reports,  and  signal-to-ncise  ratio  effects  are  estimated  to  have  been 
negligible.  The  log-normal  distribution  model  fits  the  signal  level  data  extremely  well,  at  both  frequencies; 
the  multipath  delay  spread  and  Doppler  spread  values  were  normally  distributed.  Only  the  Doppler  spread 
measurement  is  seriously  in  question.  It  is,  apparently,  tainted  by  the  60  Hz  noise  evident  on  die  frequency 
spectra. 

1.7.2  Modes  of  Propagation 

Three  modes  of  propagation  were  investigated:  turbulent  volume  scatter,  rain,  and  ducung.  Turbulent 
scatter  accounted  for  90%  of  the  hourly  median  values,  with  the  remainder  primarily  dominated  by  rain 
scatter.  Rain  on  the  path  *  as  found  to  amplify  the  Ku  band  signal,  contrary  to  expectations  based  on 
projections  matte  horn  lower- frequency  experimental  data,  but  in  agreement  with  the  theoreucal  framework 
of  die  turbulent,  thin-layer,  volume  scatter  model.  This  agrees  with  the  results  from  the  previous  two 
years.  For  neither  of  the  months  investigated,  dunng  rainy  conditions,  was  either  the  median  signal  level 
lower  than  that  of  turbulent  scatter,  or  the  percentage  of  extremely  low  reception  levels  higher.  In  general, 
the  Ku  bond  was  more  sensitive  to  rain,  reflected  in  the  greater  change  in  Doppler  frequency  shift  dunng 
rain.  Conversely,  the  C  band  was  more  sensitive  to  ducting,  a  condition  where  extended  negative  gradients 
are  coupled  into,  with  effects  similar  to  that  of  waveguides.  During  tne  few  hours  when  this  occ  ured  - 
which  was  primarily  during  the  night  hours  •  there  is  a  significant  deviation  from  the  norm,  wherefore  it  is 
considered  an  entirely  separate  phenomena  from  turbulent  scatter.  The  C  band,  and  the  Ku  band  to  a  lesser 
extent,  is  strongly  amplified  during  ducting. 

1.8  Conclusions 

This  study,  along  with  the  previous  results  of  this  same  data  collection  campaign,  show  conclusively  that 
the  Ku  band  is  reliable  for  troposcaner  communications  usage.  The  thin-layer,  turbulent  volume  scauer 
model  accurately  predicts  the  results  in  clear  air,  both  in  terms  of  total  power  received  and  the  signal  delay 
profile.  Rain  is  shown  only  to  enhance  the  signal  levels,  but  was  also  accompanied  by  a  slight  Doppler 
spread  increase.  It  is,  however,  less  severe  than  the  seasonal  Doppler  spread  vanauon.  The  wideband  system 
has  proven  very  useful  for  a  detailed  study  of  the  delay  profile;  this  aids  significantly  in  evaluaung 


6 


theoretical  models  and  investigating  weather  phenomena.  An  immediate  advantage  of  the  wideband  system 
is  that  it  allows  a  more  accurate  evaluation  of  the  categorization  of  data  into  the  various  propagauon  modes. 
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Chapter  2 

Tropospheric  Propagation  Models 

2.1  Absorption 

The  fundamental  reason  for  the  difficulties  encountered  in  the  use  of  microwave  frequencies  is  the  increased 
attenuation  experienced  in  precipitation  because  of  comparable  relative  sizes  of  the  electromagnetic 
wavelengths  and  the  water  droplets.  Figure  2.1  indicates  the  increased  attenuation  due  to  absorpuon  and 
scattering  at  the  higher  frequencies.  For  this  reason  it  was  thought  impracucal  to  use  the  Ku  band  for 
beyond-ihe-horuon  communication. 
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FIGL'RE  2.1: 


Attenuation  due  to  gaseous  consutuents  and  precipitation  [CC1K  IvQl) 
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Scattering  of  the  radiowaves  causes  attenuation,  as  does  the  absorbuon  caused  by  water  drops  with  a  ios^ 
dielectric.  This  attenuation  due  to  precipitauon  is  added  to  the  attenuation  due  to  the  gaseous  consutucnis 
in  the  air,  also  a  frequency  dependent  effect,  to  form  the  total  attenuation  loss.  The  most  accurate  estimates 
of  attenuation  by  gases  (pnmarily  water  vapour  and  oxygen)  and  by  moderate  to  heavy  rainfall  are  derived 
empirically,  taking  into  account  such  factors  as  rainrate  and  angle  to  the  horizon.  Good  theoreucal 
approximations  of  the  attenuation  have  been  reached  for  smaller  -  and  thus  more  sphencally  shaped  - 
raindrop  sizes  (Doviak  1984,  Sect.3.3].  The  attenuation  is  included  in  the  RADAR  equauon  as  integrals 
over  the  distance  of  the  paths  to  and  from  the  scanerers. 

2.2  Refractive  Index 

It  is  the  refracuve  index  (n)  fluctuations  (or,  ihose  of  the  relative  permittivity,  which  is  directly  related  to  m 
in  the  air  which  cause  the  refraction  or  reflection  of  radiowaves  back  toward  the  earth.  For  the  clear  air  these 
fluctuations  arc  caused  by  turbulent  mixing  of  the  mean  gradient  of  the  potential  refractive  index  (reflecting 
changes  due  to  the  accompanying  adiabatic  changes  in  temperature),  and  in  cases  of  precipitation  by  the 
gradient  of  the  refractive  index  across  the  different  media. 

The  complex  refractive  index  of  rain  is  expressed  as  m  *  n  *  jk,  where  k  is  the  absorpuon  coefficient,  a 
function  of  wavelength  and  temperature.  A  determination  of  the  refractive  index  is  the  first  step  toward 
finding  the  backscauer  cross  Section  of  a  particular  phenomenon. 

The  refractive  index  of  the  atmosphere  depends  pnmanly  on  the  water  vapour  pressure  which,  unlike  the 
composition  of  dry  air.  is  highly  variable: 

(1) 

N  =  77.62.  +  3.73x  103*-  :  P(mb);  T(K):  e(mb) 

T  T 

where  e  is  the  water  vapour  pressure  and  the  radio  refracuvity  N  =  (n-l)xlO6,  so  defined  because 
numerically  the  fluctuations  of  n  tend  to  be  extremely  small.  The  dependence  on  vapour  is  due  to  the 
permanent  dipole  of  the  H20  molecule,  heavily  infiuencing  the  propagation  of  short  waves.  These 
quantities  are  not  independent  of  each  other,  pressure  charges  lead  to  temperature  vanauons.  and  ihe  vapour 
pressure  is  proportional  to  the  pressure  p  by  e  =  1  62pq,  where  q  is  the  concentrauon  of  water  vapour. 
Signals  received  will  then  be  negligible  unless  temperature  perturbauons  are  very  strong  when  u,aier  vapour 
contributions  are  small.  This  dependence  on  the  water  vapour  content  explains  much  of  the  observed 
seasonal  variability  and  climate  effects  on  the  power  received. 
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Only  if  there  is  a  sufficiently  rapid  spatial  fluctuation  of  the  refractive  index  will  the  electromagnetic  energ;. 
bend  down  toward  the  earth,  where  dN/dii<-l57  per  km.  In  addition  to  the  small-scale  variations  due  to 
turbulence  there  is  a  systematic  smooth,  mean  variation  of  about  -4dN/km  (varying  from  day  to  day,  as 
well  as  seasonally)  in  the  lowest  part  of  the  troposphere,  which  in  effect  raises  the  honzon  so  that  the 
effective  earth  radius  becomes  4/3  times  the  earth's  radius,  by  straight  application  of  Snell's  law. 

The  magnitude  of  the  ratio  of  the  discontinuity  of  the  refractive  index  to  the  wavelength  can  also  determine 
the  mode  of  propagation.  For  instance,  because  higher  gradients  would  be  needed  for  short  wavelengths  10 
be  reflected,  there  is  a  dependence  on  the  ratio  such  that  reflections  are  less  likely  to  occur  in  clear  air 
condiuons  at  the  shorter  wavelengths.  It  is  therefore  possible  to  relate  the  mode  of  propagauon  to  a 
function  of  ihe  characteristics  of  the  turbulence  and  the  transmission  wavelength,  given  assumptions  about 
us  formation  and  structure.  In  this  particular  case,  at  the  higher  frequency  the  mode  of  propagauon  would 
be  turbulent  volume  scatter.  As  is  common  practice  and  most  consistent  with  present  theory,  results  were 
divided  according  to  the  modes  of  turbulent  scatter,  ducting  and  rain  scatter. 

2.3  Turbulence  in  the  Troposphere 

2.3.1  General 

By  the  late  1960's  it  was  generally  accepted  that  clear  air  turbulence  accounted  for  the  low-intensity  signal 
received  wuh  tropospheric  scatter  systems,  while  short-term  high-intensity  reception  occured  as  the  result  of 
other  mechanisms  [Tatarskii,  1971,  pl52]  such  as  ducting,  reflecuon  and  precipitauon  scaucr. 
Interpretations  of  radar  backscatter  observations  had  been  attempted  as  early  as  the  late  i930's  [Friend  1939. 
by  Ottersien.  1969].  In  the  early  60’s  explanations  of  such  "Angel  scatter"  focused  on  smooth  refracuve 
index  discontinuity  reflections  [Atlas,  1964,  by  Ouersten],  Because  of  the  diffuse  nature  :-f  the  echoes, 
how  ever,  it  was  realized  that  irregular,  turbulent  fluctuations  in  the  refracuve  index  as  a  volume  scalier 
phenomena  could  at  times  be  a  more  appropriate  model  [Atlas  and  Hardy,  1966).  [Tatarskii,  1961.  1971 
developed  mathematical  theory  of  turbulence  effects  on  radiowave  propagation.)  It  was  determined  that  a 
theory  proposing  that  the  troposphere  was  composed  of  turbulent  "blobs”  whose  spatial  extent  would  be 
described  by  an  exponential  correlation  function  [Booker  and  Gordon,  1950  p  403)  did  not  fully  explain 
experimental  RADAR  results;  the  same  experiments  seemed  to  confirm  the  correctness  of  the  turbulent 
volume  scatter  theory  for  the  high  frequency  range.  Instead  turbulence  appeared  to  be  associated  with 
stratified  layer  structures  [Richter  1969,  pl266j.  On  the  basis  of  comparisons  of  model-to-expenment 
results,  it  was  concluded  that  specular  reflections  from  these  layers  are  less  important  at  the  higher 
frequencies  where  the  turbulent  layeT  scadenng  mode  dominates,  with  transiuon  frequencies  roughly  in  the 
low  GHz  range  (Crane,  1981).  RADAR  studies  conclusively  confirmed  die  Kolmogorov  theories  ot  the 
turbulence  spectrum  upon  which  Tatarskii's  studies  and  the  idea  of  turbulent  volume  scatter  were  based 
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[KropfU,  1968).  Mulufrequency  tropo  scatter  findings  varied  greatly,  however,  and  they  in  no  way  showed 
such  conclusive  results,  calling  into  question  not  the  Kolmogorov  turbulence  energy  spectrum  or  Tatarskii  s 
theories,  but  rather  the  place  of  those  theories  in  analysis  of  trans-honzonal  phenomena.  For  this  reason 
the  idea  of  specular  reflections  remained  a  viable  altemauve  to  thin  layer  volume  scattering,  although  no 
better  at  predicung  results  at  the  higher  frequencies;  thus,  the  data  collected  at  lower  frequencies  •  by  the 
National  Bureau  of  Standards  -  is  sull  the  basis  for  models  of  the  higher  frequency  paths  used  by  the 
International  Radio  Consultative  Committee  (CCER)  [CCIR,1990,  Rice  ei.al,  TechTi'otelOl).  As  recently 
as  the  80’s,  Gossard  [1984]  has  proposed  that  layer-reflections  may  be  caused  by  stacked,  thin,  volume- 
filling  layers  with  more  complicated  and  unpredictable  frequency  dependence,  again  lending  credence  to  the 
simpler  reflecuon  model  even  at  the  higher  frequencies.  Other  directions  have  been  investigated  as  well, 
such  as  fractal  modelling  of  the  random  volume  scattering  media  [Rastrogi,  1989),  without  much  impact  m 
the  field.  Skolruk  (1974)  has  a  review  of  the  different  models.  As  many  felt  that  Tatarskii’s  ideas  were 
accurate  (because  of  agreement  with  RADAR/ Aircraft  simultaneous  measurements),  and  since  a 
mulufrequency  experiment  had  not  yet  been  performed  at  the  frequency  range  deemed  appropriate  ^previous 
higher- frequency  studies  have  focused  on  the  high  level  fields  for  the  interference  proolem),  the  present  study 
is  important  to  provide  an  indication  of  the  importance  of  the  range  of  frequencies  tested.  Previous  study  in 
this  present  data  collection  campaign  were  partly  a  further  investigation  into  these  theories  on  the  questions 
of  layer  reflecuon  vs.  turbulent  scatter.  It  was  the  finding  that  indeed  the  frequency  dependence  described  by 
the  turbulent  volume  scattering  model  is  consistent  with  the  results  [Crane  1988). 

It  is  only  with  the  complete  understanding  of  these  processes  that  a  full  description  of  the  path  loss 
predictions  of  the  RADAR  equation  (predicting  the  relation  of  the  received  to  the  transmuted  power)  can  be 
achieved.  This  theoretical  prediction  is  important,  since  empirical  predictions  for  the  higher  frequency 
ranges,  relying  on  lower  frequency  measurements,  prove  inadequate.  Remaining  important  quesuons  related 
to  turbulence  structure  include  both  the  amount  of  beamfilling  and  the  frequency  of  occurrence  of  scattering 
layers. 

2.3.2  Turbulence  Scales 

Turbulence  is  commonly  thought  of  as  the  result  of  windshearing  and  buoyancy  forces  causing  mixing  of 
air  masses  of  different  types  (i.e.with  different  temperature,  humidity  etc.).  Large-sized  inhomogeneiues 
break  up  into  smaller  and  smaller  discontinuiues,  or  eddies,  such  that  the  velocity,  temperature,  humidity, 
pressure  and  dielectric  constant  characteristics  of  such  regions  become  random  functions  of  postuon  and 


ume.  Kolmogorov  [by  Tatarskii,  1971,  pp.48-]  proposed  that  for  large  Reynold's  numbers'  the  motion 
can  be  considered  as  random,  incompressible,  flow  at  low  speeds  (compared  to  the  speed  of  sound)  and  that 
the  small-scale  velocity  field  can  be  considered  homogenous  and  isotropic.  With  continued  relative  motion 
of  the  different  air  masses  creating  eddies  of  large  spatial  scales,  there  will  be  a  spectrum  of  eddy  sizes  at  am 
one  time,  with  kinetic  energy  continually  redistributed  from  the  large  to  the  smaller  scales.  The  region  of 
the  spectra  containing  the  large  eddies  (on  the  order  of  100  m)  -  with  "outer  scale  of  turbulence '  Lo  - 
contains  the  greatest  density  of  the  average  kinetic  energy  (per  unit  mass).  The  definition  of  this  parameter 
vanes  widely.  Tatarskii's  thin  layer  outer  scale  is  a  function  of  the  turbulent  eddy  dissipation  rate,  the 
thickness  of  the  layer,  the  velocity  changes  across  the  Layer,  the  mean  wind  shear,  and  the  covariances  ot  me 
volocity  and  turbulent  heat.  The  smallest  region  (on  the  order  of  1  mm)  with  an  inner  scale  of  turbulence 
represents  the  bulk  of  energy  dissipation  beyond  which  the  destructive  action  of  viscosity  and  diffusion  have 
nullified  most  turbulent  activity: 


(2) 


4/J 

L,  =  y  ;  v  =  kinematic  viscosity  of  fluid 


Kolmogorov  proposed  in  1947  that  the  turbulence  structure  tn  the  inertia]  subrange'  well  m  between  these 
two  scales  is  dependent  only  on  the  rate  of  energy  dissipation  per  unit  mass.  £,  equalling  the  rate  of  the  flux 
of  energy  loward  Lie  smaller  scales.  Scales  in  this  range  arc  too  large  for  viscocity  effects  to  be  nouceable 
The  energy  spectrum  function  then  becomes,  by  dimensional  analysis: 


E(K)  *  Ae^k 


5  n 


;k»2fl. 
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where  A  is  a  constant  [Doviak  1981,  p346}.  Wavelengths  for  RADAR  and  satellite  communication  often 
lie  in  this  inerual  subrange,  as  do  the  vertical  wavenumbers  used  for  the  Ml  Tug -Prospect  Hill  experiment 
(See  Figure  2.2).  Note  that  for  tropospheric  scattering,  which  occurs  at  an  angle  to  the  horizontal,  the 
apparent  wavelength  in  the  z -direction,  the  factor  under  consideration  here,  is  of  a  much  greater  value  than 
the  transmission  wavelength,  but  is  still  within  the  inertial  subrange. 


Re  =  *L 

v 

where  v  and  L  are  the  charactensuc  velocity  and  scale  of  flow, 
and  v  ls  the  kinematic  vijeocity 


2.3.3  Effect  of  Turbulence  on  Wave  Propagation  and  the  Structure  Constant 

For  this  model  of  turbulence  formation,  the  propagauon  mode  and  transiuon  frequency  discussed  earlier  are 
believed  to  be  found  approximately  by  [Crane  1081.  p.656]: 

(41 

— k —  >  Lo  ;  9  is  the  scattering  angle 

2sin  0/2 

where,  if  the  inequality  is  true,  specular  reflections  are  expected  to  dominate,  and,  if  the  converse  is  true . 
thin  volume  scattering  dominates.  Due  to  the  shorter  wavelength  at  the  higher  frequencies,  the  refracve 
index  variations  span  several  Fresnel  lengths.  In  that  case  stauso.cal  correlauon  modeling  is  more 
appropriate  than  the  plane  wave  retlecuon  models  used  for  the  longer  wavelengths  -  *here  the  vanaucr.s  m 
n  are  smaller  than  the  Fresnel  length  [Crane  1988].  According  to  this  model,  volume  scauenng  should 
dominate  the  clear  air  turbulence  propagauon  on  the  Prospect  Hill  to  Ml. Tug  path  since  the  frequencies  are 
high  and  only  a  narrow  range  of  the  wavespace  (centered  around  L  *  0.5X)  is  expected  io  contribute  to  the 
scattering. 


*’  '  ^ .v» 


. . 

Figure  2.2  One -dimensional  spectrum  of  refracuve  index  fluctuations.  Note  that 

the  lower  wavenumber  for  the  ineraal  S’brange  will  move  in  an 
unpredictable  fashion  according  to  variations  in  the  turbulent  structure.  [Crane  198  i  I 

The  spaual  vanauons  of  the  fluctuations  in  the  refracuve  index  are  usually  represented  not  bv  the  special 
energy  funcuon  given  earlier  or  the  related  autocorrelation  function,  but  rather  in  lerms  of  ihe  structure 
function  defined  as: 


Dn(r^<[ri'(?i-*r>-n(?i)j2> 
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where  r  is  the  distance  between  two  measurement  points,  n'  the  deviation  from  the  mean  value  of  the  index 
of  refracaon,  and  the  symbols  <>  denote  an  ensemble  average.  With  the  earlier  assumptions  ot  isotropy 
and  homogeneity,  the  structure  funcuon  ceases  to  be  the  function  of  a  vector  distance.  This  funcuon  is  used 
since,  by  virtue  of  its  being  a  first  difference,  it  filters  out  slow  variations  in  the  fluctuation  measurements, 
as  long  as  the  distance  r  is  not  too  great.  This  fact  is  important  to  non-stauonary  processes  such  as  the 
spauai  fluctuation  of  the  refractive  index  in  the  atmosphere.  In  other  words,  it  is  usable  even  though  large- 
scale  irregularities  are  neither  homogeneous  nor  isotropic  in  the  turbulent  troposphere.  Further,  it  allows 
the  following  simple  formulation  describing  turbulence  in  the  inertial  range  (Kolmogorov  194',  by  Dew-an. 
1980,  pp.9-2 1  ] : 

(6) 


Dn(r)  -  Or*3 


This  is  Kolmogorov's  original  result,  based  on  the  assumptions  that  the  potenual  temperature  and  specific 
humidity  are  conserved  and  passive.  Cn  is  called  the  structure  constant,  defined  as: 

(7) 


In  other  words,  the  structure  constant  provides  a  measure  of  the  refractive  index  fluctuations  relauve  to  the 
eddy  sices.  It  is  directly  related  to  the  energy  dissipation  rate  by: 


C«  *  ae^3  ;  a- 1 


(8) 


where  the  energy  dissipation  consists  of  both  heating  and  diffusion  factors.  The  structure  constant  as  a 
function  of  height  thus  provides  a  basis  for  comparison  of  the  intensity  of  turbulence  as  given  by  the 
Kolmogorov  spectrum,  measured  by  different  means  with  spatially  different  data  collecuon  capabihues.  e  g 
by  RADAR.  Sondes,  airplanes. 

The  throe -dimensional  index  of  refraction  variance  spectrum  derived  from  the  Fourier  integral  ot  the  gradient 
of  the  structure  function  of  Eqn.6  over  the  enure  wave  number  space  yields: 

v9> 

$„(K)  =  0.033C„;Kun  ;  (iO<krm(rFrKr)]<23-) 

U  D 


Here: 


asdf 


K  =  k[rn(D-n(ri)  is  ihe  sexetenng  vector. 


Note  the  strong  dependence  on  the  scattering  angle.  This  important  result  is  what  enables  the  interpretation 
of  the  resolution  of  the  power  signature  at  the  different  delays  to  correspond  to  the  vertical  structure  of  the 
turbulent  layers.  In  the  one-dimensional  case,  this  dependence  on  the  wavenumber  becomes  of  ihe  power 
-5/3.3  Note  that  this  structure  applies  to  the  refractive  index  spectrum  as  a  derivative  of  the  fact  Lhat  the 
temperature  as  well  as  the  velocity  fluctuations  exhibit  this  spectral  wavenumber  dependence. 

This  result  can  then  be  applied  toward  the  formulation  of  the  scattering  cross  Secuon  for  a  final  defimuon 
of  the  RADAR  equation.  From  general  Bragg  scatter,  the  resulting  wavelength  dependence  of  the  cress 
Secuon  is  [Doviak  1984,  p385]: 

t.10) 

q  =  87iVOn(K)m  ;m  =  sin2*  is  a  polarization  mismatch  factor 

;X  is  ihe  angle  of  the  polaruauon  of  the  incident 
wave  to  the  scattering  direcuon 

Because  its  value  is  near  unity  at  the  angles  under  c  > .  iderauon,  the  polarization  facior  is  ignored  in  the 
computer  models  for  the  path.  These  formulas,  along  with  correct  modelling  of  the  propagation  path 
geometries  and  accurate  assumptions  about  the  general  structure  of  turbulence,  establishes  the  groundwork 
for  the  relation  of  the  reception  characteristics  to  the  intensity  of  clear  air  turbulence. 

In  the  final  version  of  the  RADAR  equation  this  wavelength  dependence  is  of  the  power  - 1/3.  Previous 
experiments  have  yielded  ambivalent  results  as  to  this  frequency  dependency,  ranging  from  - 1  to  3  fEklund 
and  Wickerts.  1968],  As  Figure  2.2  shows,  however,  it  is  important  to  ensure  that  the  frequencies  lie  well 
within  the  inertial  subrange  of  the  Kolmogorov  spectrum  for  the  turbulent  volume  scatter  to  apply.  The 
boundary  between  the  buoyancy  and  inertial  subranges  is  not  a  fixed  one:  rather,  it  is  dependent  on  the 
particular  weather  conditions.  This  project  has  represented  the  first  such  study  since  ihe  many  previous 


m(r)  *  unit  vector  from  the  tiansmitter  to  ihe  scatterer. 
n(F)  =*  unit  vector  from  the  scalier er  to  ihe  receiver, 
such  that  they  define  a  scattering  angle  making: 

K  =  2ksm&-  .  where  again  k  =  *3- 

2  a. 

The  structure  function  has  been  relaxed  to  ihe  eneTgy  of  the  t  rbulence.  and  ihe  spectrum  of  the 
structure  funcuon.  (V  m  the  one-dimensional  case)  is  related  to  die  spectral  energy  density  E  by: 

27tk2  2nk  dk 


studies  nave  been  in  too  low  a  range,  the  frequencies  lying  sow  in  the  inertial  subrange,  nov  m  the 
buoyancy  subrange. 

2.3.4  Turbulence  Structure 

Several  studies,  by  use  of  radio-sondes  and  aircraft  as  well  as  RADAR,  have  clarified  the  strucnire  of 
turbulence  in  the  troposphere.  They  show  that  the  enure  troposphere  can  be  thought  of  as  containing  a 
weak  background'  turbulence.  High-intensity  layers,  measuring  from  tens  of  meters  to  a  few  hundred 
meters  thick.  a..d  of  up  to  several  tens  of  kilometers  in  horizontal  extent,  lie  within  this  "background 
turbulence.  These  high-intensity  layers  are  much  more  likely  than  the  background  turbulence  to  be  detected 
by  RADAR.  They  are  less  often  detected  by  aircraft,  however,  and  sondes  never  reveal  their  honsontal 
extent.  These  layers  contribute  a  large  poruon  of  the  scattered  signal  in  troposcaaer  systems  during  clear  ar 
condiuons,  depend  ,.g  on  the  vertical  extent  of  the  layers,  their  spatial  placement,  and  the  strength  of  the 
turbulence-generating  shearing  force*.  Appearing  randomly  over  ume,  in  uniform  spaual  distnbuuon,  the 
layera  cause  every  part  of  the  c,  a  volume  to  be  turbulent  at  some  ume  (Crane  1980).  For  this  reason 
-  and  since  any  one  wavelength  is  responsive  only  to  a  particular,  narrow,  range  of  scale  of  turbulence  -  the 
strong  layering  effect  is  expected  to  occur  only  within  a  limited  volume  a  small  percentage  of  the  ume. 

The  occurrence  of  particular  structures  should  be  taken  into  account  for  proper  prediction  of  signal  recepuon 
in  clear  air.  For  instance,  RADAR  echoes  often  show  wave  patterns  in  clear  air.  In  the  common  case  of 
vertical  wind  shear,  such  waves  would  presumably  exhibit  stronger  refracuvity  variance  in  us  upper  and 
lowei  extremes  when  breaking  waves  are  formed  (Ottersten  1969).  This  may  only  be  detected  as  two 
distinct  layers  to  a  troposcaaer  system.  Along  with  ducting,  these  cases  are  more  pertinent  to  the 
understanding  of  interference  phenomena.  Thermals,  or  convecuve  bubbles,  are  expected  to  reach  the  lower 
edge  of  the  common  volume  at  about  ten  am. 

2.4  Ducting 

Although  ducting  is  a  special  case  of  clear  turbulence,  it  is  usually  not  considered  as  :uch  During  ducting, 
layers  with  particularly  high  dn/dh  boundaries  form  parallel  planes;  sometimes  microwaves  can  enter  the 
guides  and  become  trapped,  in  a  waveguide-hke  fashion.  A  much  larger  portion  of  the  transmitted  signal 
will  w  this  case  be  directed  toward  the  receiver.  Because  of  the  unsteady  nature  of  the  ducts  and  the  fact  ’ha; 
they  will  probably  not  be  quite  planar  makes  modelling  of  this  phenomenon  very  difficult.  curther 
complexity  is  introduced  since  the  coupling  into  the  duct  of  the  different  modes  -  and  the  resulting  spectral 
dependence  -  is  difficult  to  predict.  In  genera!  the  lower  frequency  is.  as  one  woo'd  expect,  considerably 
more  sensitive  to  this  high-level  power  phenomenon.  (Further  reducuon  of  the  frequency,  however 
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the  verucal  extent  of  the  duct  to  become  more  important.)  From  observations  of  long  ume-lags  between 
the  constructive  and  destructive  interference  patterns  that  sometimes  exist  in  the  waveguide  [Crane  198  i ). 
the  ducting  mean  Doppler  is  expected  to  be  very  low. 


2.5  Precipitation 


There  is  more  detailed  knowledge  of  rain  formation  and  characteristics:  its  drop  size,  shape  and  distnbuuon 
effects  on  radio'  ve  propagation  are  comparatively  well  understood,  primarily  because  of  Lhe  relative  ease  in 
comparing  data  collected  by  different  means  and  the  simpler  descriptive  nature  of  this  phenomenon.  Ram 
formation  during  the  different  periods  of  the  year  has  also  been  studied  extensively,  and  snow  and  ice  panicle 
distribution  can  be  predicted  with  some  accuracy  at  different  heights  from  ram  measurements  and  height 
profiles  of  temperature  and  humidity.  Rain  is  eminently  detectable;  it  is  easily  differentiated  from  other 
forms  of  scatter  by  a  large  positive  Doppler  signature.  (We  have  defined  posiuve  Doppler  as  that  which 
shortens  the  wave  paihiengih.  Rain,  with  a  downward  motion  toward  the  straight  line  passing  through!  the 
earth  that  connects  the  two  antennas,  has  precisely  this  effect.)  The  backscatter  cross  Section  of  a  water 
drop  per  unit  volume  •  Lhe  Rayleigh  approximation  -  is  (Crane  1974]: 


;D-Drop  diameter 
;K*  » (rnM)/(m2*2) 


01) 


where  m  is  the  complex  refractive  index  of  water,  and  IKwl  is  about  0.96  for  water,  0.42  for  ice  spheres. 
The  total  scatter  from  rain  and  water  droplets  will  then  simply  depend  on  the  refiecuvity  factor: 

02) 

NoD‘dD 

, 

where  N  is  the  number  of  drops  in  a  unit  volume  This  wtU  be  modified  later  by  other  considerations 
Now  an  important  di<ference  in  the  scattering  cross  Sections  of  clear  air  turbulence  and  ram  scatter  become  s 
evident;  while  turbulence  has  a  very  strong  dependence  on  the  scattering  angle,  ram  scatter  is  independent  .-i 
the  angle  for  linear  polarization  perpendicular  to  the  plane  ot  scattering,  or  propc  onal  to 
cos^ft)  for  lincai-  polarization  m  the  plane  oi  scattering. 


The  uoposcaucr  system  is  clearly  more  sensitive  to  tain  scalier.  The  greater  attenuation  experienced  m 
precipiiauon,  however,  causes  the  received  power  in  nun,  relative  to  that  ji  clear-air  conditions,  to  depend 
on  the  intensity  of  the  rain.  This  trade-off  is  very  important,  it  was  expected  that  ram  on  the  path  wou!J 
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simply  attenuate  the  signal  so  much  that  it  would  decrease  the  reliability  of  the  link.  It  has  been  shown, 
however,  that  the  trade-off  is  sufficient,  so  that  statistically  the  link  is  no  less  reliable  during  precipitauon. 


A  comparison  of  the  scatter  cross  Sections  of  rain  and  turbulent  clear  air  condiuons  shows  (Figure  2.3)  that 
extremely  clear  air  conditions  are  required  for  turbulence  to  provide  a  significant  portion  of  the  scattering 
cross  Section  and  thereby  dominate  the  mode  of  scattering.  Because  this  plot  is  for  RADAR  at  a  scattering 
angle  of  180°.  compensation  should  be  made  for  cropo  scat  ter  for  small  angles  due  to  the  effectively  enlarged 
scattering  cross  Section  of  turbulence.  Therefore,  the  proper  placement  of  the  lines  is  about  20  dB  higher 
than  in  the  Figure  shown. 
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Figure  2.2:  Plot  of  the  structure  constant  under  conditions  for  which  scatter 

contributions  from  CAT  and  rain  are  expected  to  be  equal  (Go»sard  19V 1 1. 
Note  that  there  is  a  required  upward  shift  of  20  dB  of  ihe  lines  for  the 
Prospect  Hill-Mt-Tug  path. 
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2.6  Terrain  Diffraction 

Terrain  diffraction  plays  an  important  role  in  tianshorizon  communication;  it  cun  even  dominate  me 
propagation  process.  Diffraction  is  most  important  when  there  is  a  common  edge  (such  as  a  mountain 
peak)  at  the  horizon,  located  at  a  very  small  angle  to  both  the  transmitter  and  receiver.  Under  the  right 
conditions,  it  is  possible  that  diffraction  from  several  consecutive  edges  can  dominate.  Further,  during 
times  with  unusually  high  refractive  index  variation,  the  radio  horizon  can  be  altered  to  a  high  degree  so  as 
to  lower  the  angles  with  any  distant  edges.  The  Prospect  Hill  -  Mt.Tug  path  is  not  expected  to  be 
significantly  affected  by  this  phenomena.  The  angles  with  the  honzon  are  large  (see  footnote  5  for  a  lisung 
of  some  important  parameters),  and  even  during  super-refractivity  conditions,  the  waves  must  diffract 
around  several  edges.  If  diffracuon  occurs,  there  would  be  an  increase  of  the  power  at  the  lowest  height 
viewed  by  the  system  (corresponding  to  the  shortest  path).  The  Doppler  shift  associated  with  that  power 
would  also  be  very  small.  Neither  of  these  characteristics  are  observed  for  any  extended  penod  of  ume  on 
the  path. 


2.7  Aircraft  Scatter 

The  most  important  anomalous  form  of  scatter  is  from  aircraft  Aircraft  in  the  common  volume  on  the 
Logan-Lebanon  route  are  detected  clearly  by  rapidly  vanishing  high-level  fields  along  with  a  strong  Doppler 
signature.  By  analyzing  the  data  in  terms  of  median,  rather  than  mean,  quantities,  such  shot  noise  is 
expected  to  have  minimal  effect 
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Chapter  3 

The  System  and  its  Channel  Characteristics 
3.1  General 

Troposcaner  systems  consist  of  a  transmitter,  the  tropospheric  scatter  channel,  and  a  receiver  located  bc>ond 
the  radio  horizon.  A  high-gain  receiver  antenna  picks  up  the  minute  amount  of  transmitted  energy  scattered 
towards  it  from  the  atmosphere.  The  directions  in  which  the  antennas  point  most  often  define  a  great  circle 
path  (the  shortest  distance  path).  Both  dish  antennas  are  highly  directive,  such  that  a  small  common 
volume'  is  defined  within  their  circular  3-dB  beamwidihs.  Because  of  the  low  pointing  ot  the  antennas  th  s 
common  volume  is  roughly  cylindrical  along  the  great  circle  path.  In  actuality,  although  the  beam  leaving 
the  transmitting  antenna  is  of  an  in-phase  cylindrical  shape  mimicking  the  antenna  face,  the  beam  >preaJi  to 
a  slightly  conical  shape  as  it  is  diffracted  in  the  atmosphere.  For  our  particular  path  the  Ku-frequency  i ! ;  c. 
18  GHz,  though  nowadays  said  to  lie  in  the  J  frequency  range)  volume  ranged  from  1.1  km  in  height  abo^e 
sea  level,  or  400  m  above  the  mountain  range  below,  to  2.3  km  above  sea  level,  and  horizontally  covered 
roughly  60  km,  over  a  161  km  long  path  between  Mt.Tug  in  Lebanon,  NH,  and  Prospect  Hill.  Waltham, 
MA4.  bince  the  common  volume  is  closer  to  the  Mt.Tug  side  and  the  antenna  is  for  the  Ku-band  of  the 
same  size,  this  antenna  defines  the  limits  of  the  common  volume.  Simultaneous  measurements  were  made 
at  a  lower,  more  commonly  used  frequency  in  the  C-band  (4-8  GHz.  nowadays  in  the  G-band),  with  larger 
antennas  resulting  in  a  smaller  common  volume  ranging  from  1.1  to  1.7  km  in  height  and  40  km 
horizontally.  (The  beam  width  is  given  by  0  =  aX/D,  where  a  is  a  constant  determined  from  considerauon, 
of  antenna  shape  and  pattern.)  The  axes  of  the  antennas  were  pointed  in  the  same  directions  so  that  the 
experiment  is  a  multi-frequency  measurement  of  roughly  the  same  volume.  Not  only  the  scattered  within 
the  volumes  defined  by  the  half-power  points  are  important;  energy  retransmitted  from  die  enure  volume 
encapsulating  these  two  cransmit-receive  points  and  beyond  also  reaches  the  receiver,  albeit  weighted  by  the 
antenna  pattern.  Of  course,  it  is  also  affected  by  the  aforementioned  relations  of  the  different  scattering 
modes. 


Prospect  Hill  -  Mt.Tug  Bis  tine  Scatter  System 


Frequency 

15.731  GHz 

4  95  GHz 

Distance 

160.9  km 

Trens.  Height 

159.4  m 

Rec.  Height 

493.8  m 

Trans.  Aperture 

3  ft 

29  ft 

Trans.  BW 

26.49  mrad 

8.36  mrad 

Rec  Aperture 

3  ft 

10  ft 

Rec.  BW 

26.49  mrad 

25.25  mrad 

Scattering  Angle 

35.23  mrad 

34.61  mrad 

Modulation 

400  MBit/Sec 

None 

Carrier  Power 

225  W 

30  W 

Rec.  Noise  Level 

•144  dBm 

-147  dBm 

Several  factors  other  than  the  types  of  scattering  described  earlier  affect  the  magnitude  of  the  received  power. 
These  include  the  effects  of  the  antenna  pattern  and  the  so<alled  'aperture-to-medium'  coupling  as  well  as 
the  path  geometries. 


The  far-field  antenna  power  density  patterns  art  derived  from  the  two-dimensional  Founer  transforms  of  the 
field  in  the  antenna  aperture,  for  the  case  of  the  paraboloid  reflectors  resulting  in  (normalized,  Gt  will  denote 
the  magnitude  along  the  axis): 


f2(0)  _  |8JMtDsin9)A.]|* 
1  [(TtDsin6)A]2  1 


;  $=angular  distance  from  beam  axis 
;  D=antenna  diameter 
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The  antenna  pactem  then  consists  of  a  main  beam  and  a  number  of  sidelobes  well  outside  of  ihe  common 
volume.  As  a  result,  there  is  severe  suppression  of  the  power  received  from  scaoerers  outside  the  common 
volume,  though  this  suppression  is  not  uniform.  Nevertheless,  it  ts  conceivable  that  strong  enough 
scaoerers  outside  the  common  volume  relative  to  scaoerers  within  the  common  volume  could  significantly 
add  to  or  even  drown  the  effect  of  lesser  scaoerers  inside  the  common  volume. 


Not  only  the  strength  of  the  scaoerers  but  also  the  extent  to  which  the  scaoerers  fill  ihe  volume  plays  an 
important  role.  A  large  volume  of  scaoerers  will  retransmit  more  energy  toward  the  receiver  than  a  lesser 
volume  of  otherwise  equal  scaoerers.The  far-field  antenna  gam  function  is  given  by: 

(14) 

0(8.0)  ;Ac«kAf1(e)=.65Af*(e>  ;A=aCl 

\2  * 

G(0.0)S  62^(0) 

X2 


As  antennas  increase  in  size  and  thus  decrease  in  beamwidth  to  achieve  higher  gams,  there  is  therefore  an 
accompan>  :rg  reduction  of  the  volume  of  scaoerers.  This  so-called  aperture  to  medium  coupling  loss  is 
accounted  for  by  the  standard  references  used  for  design  (CCIR,  1990.  Rice  ei.al,  Tech. Note  10 1  j  by  a  simple 
formula  relating  the  loss  to  the  sum  of  the  two  antenna  gains.  This  method  lacks,  however,  in  that  n  tails 
to  address  the  question  of  beamfiUing  or  geometric  structure  of  turbulence  or  other  scaoerers.  For  instance, 
if  a  is  found  that,  for  a  significant  portion  of  the  lime,  thin  turbulent  layers  contained  within  a  small 
poruon  of  the  common  volume  contribute  to  the  received  power,  then  it  is  reasonable  to  assume  that  use  ot 
an  aiuenna  with  higher  gain  would  not  entail  a  coupling  loss,  since  any  reduction  in  the  size  of  the 
common  volume  would  sail  contain  the  scaoerers.  A  correction  factor  taking  into  account  layer  thickness 


:i 


is  introduced  in  the  models  (Ottersten  1969,  Crane  1980,  p  185] ,  but  encounters  difficulties  if  the  structure 
is  more  complex. 

The  coupling  effect  also  plays  an  important  role  as  to  the  relative  importance  of  the  scatiercrs  inside  and 
outside  of  the  common  volume.  It  might  appear  that  for  precipitation  the  antenna  function  and  for 
turbulence  the  added  angle  dependence  given  in  equation  9  might  eliminate  all  impact  of  scattercrs  in  a 
region  outside  the  common  volume.  By  virtue  of  its  large  volume,  though,  the  summed  effects  from  such 
a  region  may  in  fact  be  greater  than  from  the  common  volume  for  the  precipitation  case.  For  instance,  if 
there  is  rainfall  direcdy  above  ihe  transmitter,  such  as  would  enhance  forward  scatter  at  these  frequencies,  a 
very  large  angular  region  of  the  transmitting  antenna  beam  would  then  be  filled  by  scaoerers.  corresponding 
to  a  large  volume  covered  by  the  receiver.  The  effect  from  the  rain,  even  when  absent  from  the  common 
volume  could  thereby  produce  a  signal  even  stronger  than  the  norm. 

3.2  Channel  Characteristics 

Along  with  the  attenuation  caused  by  absorption  and  scattering,  the  channel  causes  additional,  time-variant, 
fading  of  the  signal  level.  The  fading  is  categorized  into  short-term  and  long-term  fades.  The  rapid  fading 
(on  the  order  of  seconds)  is  the  result  of  random  phase  additions  of  the  scauerers  causing  constructive  and 
destructive  interference.  In  accordance  with  theory,  the  power  level  is  Rayleigh  distributed  with  respect  to 
the  fast  fade;  since  the  large  number  of  scaoerers  with  uniformly  distributed  phases  result  in  a  gaits., un 
process,  the  in-phase  and  quadrature  components  are  normally  distributed.  The  slow  fade  (on  the  order  of 
tens  of  minutes),  which  is  the  result  of  changes  in  the  weather,  has  been  empirically  determined  to  bo  log- 
normally  distributed.  In  terms  of  the  diurnal  and  seasonal  changes  there  is,  generally,  an  increase  in  the 
signal  level  dunng  the  warmer  hours  of  the  day,  and  during  the  wanner  and  more  humid  seasons.  The 
design  of  troposcatter  communications  systems  is  based  on  estimates  for  channel  condihons  of  the  worst 
month'  of  the  year,  new  statistical  methods  for  doing  so  have  recently  been  developed  (Crane  1990], 

Complete  characterization  of  the  channel  requires  not  only  knowledge  of  the  signal  level  vanaoon;  ihe 
Doppler  signature  and  the  delay  structure  also  provide  crucial  design  parameters.  Only  when  this 
characterization  c.  the  channel  has  been  made  complete  can  determination  of  the  bandwidth,  proper  coding, 
and  diversity  schemes  be  made. 

3.2.1  Delay  Spread 

The  pathlengih  associated  with  eroh  scauerer  results  in  a  particular  time  delay  between  ihe  ume  of 
transmission  and  reception.  The  sum  of  the  power  received  from  all  equidelay  paths  consututes  the  total 


power  received  ai  any  one  delay.  Equidelay  paths,  usually  expressed  in  nanoseconds,  form  near  ellipsoids 
around  ihe  sites.  Due  to  the  angular  dependence  -  from  the  antenna  gain  function  and  scattering  cross 
Section  definitions  -  and  a  greater  amount  of  absorption,  the  contribution  from  scaoerers  at  greater  delays  is 
sharply  diminished.  The  exact  shape  of  the  powcr-vs-delay  signature  -  the  multi-path  intensity  profile  - 
depends  on  the  spatial  distribution  and  nature  of  the  scattertrs.  Figure  3.1  a-b  shows  two  typical  clear  air 
profiles. 

A  quantitative  measure  of  the  delay  signature,  the  delay  spread  of  the  signal,  is  defined  as  ihe  distance  tin 
nanoseconds)  between  the  half-power  points  of  the  multipath  intensity  profile.  This  is  a  very  important 
quantity  since  us  reciprocal  is  a  measure  of  the  coherence  bandwidth  of  the  channel;  the  ume  dispersion  of 
the  signal  indicates  the  potential  intersymbol  interference  (IS I)  on  very  rapid  data  screams. 

For  turbulent  scatter,  the  mam  source  of  reradiated  energy  will  derive  from  the  common  volume.  It  is 
possible,  therefore,  to  consider  ihe  delays  to  represent  different  heights  in  that  region.  This  is  the  direct 
result  of  the  angular  dependence  of  the  scattered  power  given  in  equation  9.  Turbulence-scattered  power 
received  at  any  one  delay  is  the  sum  of  powers  reradiaied  from  all  points  on  a  corresponding  sphenodal 
shell;  the  planes  are  defined  by  the  intersection  of  ellipsoids  (associated  with  particular  delays)  and  the 
common  volume.  Within  the  common  volume,  horizontal  deviation  from  the  center  along  the  ellipsoid 
causes  little  weighting  by  the  relation  given  in  eqn.  9.  Vertical  deviation,  however,  leads  to  significant 
change  in  the  scattering  angle,  even  within  the  common  volume  .  For  a  uniformly  turbulent-scaaerer  filled 
volume,  computer  models  show,  the  received  power  theoretically  yields  a  vertical  profile  roughly  Uke  that 
of  a  rapidly  rising  exponential,  followed  by  a  slower  exponential  decay,  see  Figures  3  2a)  and  b).  Because 
of  this  general  shape,  the  delay  spread  is  expected  to  be  small  dunng  clear  air  conditions.  During  dueling, 
or  conditions  where  thin,  high-intensity,  layers  dominate  the  scattering  process,  the  spread  will  be  minimal. 

The  multi-path  intensity  profile  aids  also  in  discerning  the  mode  of  operation  of  the  channel.  Clearly,  if 
there  is  only  one  layer  causing  specular  reflection,  the  delay  profile  would  show  a  narrow  spike, 
characterized  by  very  slow  motion,  randomly  positioned  in  the  common  volume.  No  such  spikes  exist  if. 
however,  the  volume  consists  of  many,  thinly  layered,  reflecting  surfaces,  it  could  be  very  difficult  to 
distinguish  them  from  a  scattering  process  caused  by  turbulent  motion  by  simply  looking  at  the  multi-paih 
intensity  profile. 

It  is  not  possible  with  certainty  to  make  any  height-to-deLay  association  in  the  case  of  rain,  since  ihe 
angular  dependence  of  Equation  9  is  Lacking  in  the  cross  Section  equation  for  rain  (Equauon  11).  !i  is  at 
limes  possible,  of  course,  that  this  association  is  in  fact  accurate.  It  is  evident  from  Figure  2.3  dial  ver. 
clear  conditions  are  necessary  for  a  completely  accurate  vertical  posiuonmg  of  liie  scailercrs  by  inference 


23 


from  the  delay  profile.  The  delay  spread,  therefore,  is  very  large  during  precipitaaon.  In  fact,  under  the 
nght  conditions,  the  delay  signature  can  be  virtually  flat,  even  when  the  received  power  is  very  high. 
Computer  modelling  shows  that  a  locally  limited  volume  of  rain,  positioned  at  different  Sections  across  the 
region,  results  in  a  variation  of  received  power  of  less  than  10  dB.  There  are  some  indications  by  which  it 
could  sull  be  possible  to  estimate  where  the  rain  is  located,  if  it  is  locally  limited  at  all.  For  instance,  rain 
directly  above  the  transmitter  would  be  accompanied  by  substantial  increase  in  the  variance. 


a)  b) 

Figure  3.1:  Power  vs.  delay  signatures  during  a)  dear  air,  uniformly  filled  common 

volume,  b)  clear  air,  layering  structure. 


In  order  to  find  the  power  associated  with  each  delay  it  is  necessary  to  fust  establish  a  rclauve  ume 
reference.  To  aid  in  this,  a  109.3  bit  long  phase-modulated  pseudo-random  no>se  code  generator  was  applied 
to  the  Ku  band  transmission.  An  identical  device  ,  switching  the  phase  of  a  correlating  signal  at  the  same 
rate,  was  applied  at  the  receiver  end  The  power  received,  when  the  transmitted  signal  is  at  a  delay  such  that 
the  incoming  and  correlating  signals  are  in  phase,  is  1023  umes  the  power  at  an  out-of-phase  delay.  The 
correlating  device  is  then  clocked  to  skip  one  bit  to  form  the  next  1023  bit  long  posiuon,  and  the  power  at 
another  delay  can  be  recorded.  The  device  is  called  a  RAKE  correlator.  The  envelope  of  the  power  as  a 
function  of  delay  yields  the  multipath  intensity  profile.  Since  the  rate  of  the  PRN  code  is  JOO  MBiuscc. 
the  ume-resoiuuon  of  the  device  is  2.5  ns,  yielding  in  the  dear-au  case  a  spatial  resolution  of  roughly  30 
m.  This  is  a  much  greater  resolution  of  low  tropospheric  RADAR  data  than  has  been  achieved  before. 
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Resolution  volumes  of  similar  size  have  been  achieved  using  RADAR  [e  g.  Woodman  1980]  but  not  at 
such  low  altitudes. 


Multipath  Delay  Profile 

Calculated  tor  the  Prospect  Hitt  to  Mt  Tug  i‘r^poscatter  Path 
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a)  Theoretical  plot  of  delay  signature  of  the  signal  level,  showing  also  the  effect  of 
two  strongly  turbulent,  100m  thick,  layers,  b)  Histogram  of  the  delay  signature 
summed  over  a  week  s  observauons,  in  terms  of  number  of  meauremenis  above  the 
3-dB  points,  implying  a  uniform  distribution  of  the  thin  layers  with  height. 


Figure  3.2: 
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While  the  delay  spread  limits  the  uncoded  bandwidth  of  the  channel  as  l/(delay  spread),  coding  and  various 
demodulation  schemes  can  improve  this  rauo,  subject  to  the  limits  of  the  signal-to-noise  ratio.  The  RAKE 
correlator  is  in  fact  the  optimum  receiver  for  processing  the  wide-band  signal  (i.e.  where  bandwidth  >> 
inf..-  i  rate).  Without  coding  a  typical  delay  spread  of  150  ns  would  allow  for  a  bandwidth  of  about  7 
MHz;  with  the  RAKE  receiver  the  theoretically  achievable  bandwidth  equals  the  bit  resolution,  i.e.  the  code 
rate  or  400  MBit/sec. 

3.2.2  The  RADAR  Equation 

Power  reception  is  straight-forward  for  the  volume  forward  scatter  model,  once  ihe  scattering  cross  Section 
has  been  determined.  The  power  density  incident  on  the  scatierer  is: 

S,  *  ;x  »  position  of  scatterer  relative  to  transmitter 

4rart2 

The  power  radiated  by  the  scatterer  toward  the  receiver  is: 

S,  a  ■♦).  ;p  *  position  of  receiver  relative  to  scaoerer 
47trp2 

The  total  power  at  the  receiver  is  therefore: 

Pr  ■  SrA«  »  Gr(9,0) 

4re 

Taking  into  account  line  losses  and  the  specific  attenuauon  losses  from  the  atmosphere  the  final  equation 
then  becomes: 

05) 

P‘  (4*r  Jvd  rtfp 

where  d  is  the  location  of  the  receiver  relative  to  the  transmitter  [Crane  1981].  The  scattering  cross 
Sections  are  as  given  in  chapter  2.  This  equauon  can  be  fitted  to  yield  the  power  received  at  any  delay,  'he 
total  power  etc,  by  simply  changing  the  values  (i.e.  Cn  for  clear  air)  in  the  scattering  cross  Section.  Such 


changes  effectively  simulate  litfereni  atmosphenc  condiuons.  By  doing  so,  layered  turbulence  structures  on 
top  of  weak  background  turbulence  have  been  modeled,  as  have  different  rain  configurations. 

3.2. 3  Mean  Doppler  Frequency  of  the  Doppler  Spread 

The  Doppler  frequency  shift  in  the  received  signal  was  measured  at  each  delay  as  well  as  the  power.  It  is 
estimated  from  the  in-phase  and  quadrature  components  of  the  received  signal  by  knowing  the  frequency 
which  caused  the  shift  The  in-phase  and  quadrature  components  of  the  signal  are  created  for  this  express 
purpose.  The  Doppler  roughly  measures  the  component  of  motion  of  the  scatrerers  along  the  verucal  axis; 
it  is  estimated  by  the  pulse-pair  algorithm.  The  mean  Dopp'er  is  simply  the  power-weighted  average  of  the 
Doppler  shift  of  the  enure  volume  of  scatxerers.  Any  mouon  of  a  scatterer  -  other  than  along  the  line  ot  ihc 
great  circle  path  (or,  rather,  other  than  along  the  equidelay  planes)  -  will  in  effect  either  leng’hen  or  shorten 
the  path  of  the  scatterered  energy.  The  motion  will  cause  a  negative  or  a  positive  Doppler  shift, 
respectively,  of  the  wave.  For  a  symmetric  path,  the  horizontal  components  of  the  scatterers  will  cancel 
each  other  out;  the  motion  away  from  and  toward  the  center  cany  equal  weighting  from  the  antenna 
functions  and  the  scattering  cross  Sections,  and  in  equal  amounts  lengthen  and  shorten  the  path  of  the  wave 
(This  assumes  that  the  horizontal  wind  is  uniform  over  the  common  volume,  a  reasonable  assumption  ) 
Vertical  motion,  on  the  other  hand,  does  not  cancel  out.  Downward  motion,  for  instance,  always  has  a 
component  that  shortens  the  path,  none  that  lengthen  it.  A  clearly  positive  Doppler  is  therefore  detected 
during  precipitauon. 

Small  amounts  of  rain  may  not  have  sufficient  downward  motion  over  the  enure  volume  to  cause  a  large 
Doppler  shift.  Therefore  it  is  theoretically  possible  to  have  rain  on  the  path,  but  for  the  data  to  be 
considered  clear  air  data.  During  clear  air,  and  possibly  during  some  cases  of  precipitauon,  the  Doppler  at 
different  delays  represents  the  motion  of  the  scanerers  at  different  heights.  It  i. s  posss'o'e  to  observe  the 
process  of  various  meteorological  phenomena  using  these  signatures,  such  as  slowly  tailing  snow  melung 
to  form  rapidly  falling  rain. 

The  Doppler  spread  is  a  measure  of  the  rapidity  of  fading  on  the  path.  It  is  defined  as  the  instance  between 
the  half-power  points  on  the  power  spectrum  Like  the  delay  spread,  uiis  is  an  important  measure.  It  is 
crucial  in  determining  the  «ymg  distances  for  frequency  and  phase  modulation  schemes.  Since  the  delay 
spread  is  much  smaller  than  the  coherence  ume  (l/(Doppier  spread)),  the  fading  is  always  slow  compared  10 
the  signal  rate.  The  signal  transmission  will  therefore  be  frequency-rionselective.  It  is  important  to 
understand  that  the  Doppler  spread  is  the  frequency -spread  at  any  one  delay  (usually  the  one  associated  with 
die  greatest  power),  and  that  thr  Doppler  shift  is  the  shift  in  ihe  spectrum  at  any  one  delay.  This  o  m 
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contrast  to  the  mean  Doppler,  which  is  an  average  measure  over  the  enure  volume.  Figure  3.3  shows  a 
typical  spectrum 
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Figure  3.3:  a)  frequency  spectra  at  the  maximum  power  delay  tap,  showing  Doppler 

shift  and  spread. 

There  are  several  factors  that  together  cause  the  Doppler  spread  There  is  a  component  from  equation  3, 
since  the  total  velocity  variance  is  an  integral  of  the  spectrum  EQc).  The  Doppler  spread  would  be  directly 
proportional  'jo  the  turbulent  dissipation  rate,  if  this  was  the  only  component  Unfortunately,  however, 
there  is  additional  broadening  of  the  spectrum:  by  the  volume  averaging  of  the  scauerer  contributions  within 
any  one  delay  shell,  by  the  time  averaging  of  several  measurements,  and  by  the  wind  shear.  While  it  may 
be  possible  to  separate  out  the  variance  due  to  turbulent  activity  in  some  RADAR  measurements,  this  is 
not  possible  on  our  path  since  the  wind  shear  component  is  by  far  the  dominating  mode.  Even  m  uniform 
wind  with  no  shear,  there  would  be  a  range  of  wind  components  in  the  direction  toward  and  away  from  Lhe 
center  of  the  path. 

3.3  Particulars  of  the  Prospect  Hill  *  Mt.Tug  Pi.: 

There  are  practical  problems  associated  with  installing  a  troposcatter  system.  These  nave  resulted  in 
alterations  of  the  channel  characteristics  described  above,  unique  to  the  MLTug  -  Prospect  Hill  path.  Most 
significant  among  these  is  the  offset,  or  misalignment,  of  the  antennas.  It  is  the  result  of  the  practical 
'definition'  of  the  great  circle  path  when  using  less  sophisticated  equipment. 
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In  practice,  aligning  the  antennas  is  done  by  scanning  the  honzon  with  the  receiver  antenna  until  the 
position  corresponding  to  the  highest  power  is  found  This  is  a  time-average  operation,  and  is  subject  to 
change  with  different  weather  conditions.  Furthermore,  in  the  mountainous  terrain  between  Boston  and 
Lebanon,  irregularities  in  the  local  horizons  of  the  antennas  will  offset  even  the  true  maximum  power 
position.  The  lower  picture  in  Figure  3.4  shows  how  the  antennas  could  most  closely  match'  if  the 
receiver  was  pointed  at  a  valley,  instead  of  along  the  great  circle  path  -  where  it  would  point  at  an 
ntervemng  mountain. 

This  misalignment  or  irregularity  in  the  local  honzon  plays  an  important  role  in  the  interpretauon  of  trie 
data.  As  shown  in  the  Figure,  even  if  there  is  only  a  uniform  horizontal  wind,  the  components  of  the  wind 
will  no  longer  cancel  oul  As  shown,  westerly  and  easterly  winds  yield  negauve  and  postuve  Doppler 
signatures,  respectively.  The  upper  picture  in  Figure  3.4  shows  that,  even  if  the  antennas  are  aligned, 
intervening  mountains  can  cause  such  Doppler  shifts. 
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Figure  3,4: 


The  local  honzon  and/or  alignment  offset  of  the  antennae,  causing 
a  negauve  Doppler  signature  during  westerly  winds. 
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Figure  3.5: 
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March  tropcscauer  mean  Doppler  shift  measurements  of  lower  troposphere  and 
Concord  ground  wind  measurements 
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The  fact  that  the  component  in  the  vertical  direction  due  to  the  horizontal  wind  is  small,  is  counteracted  by 
us  greater  absolute  value;  vertical  movement  is  usually  on  the  order  of  cm/sec,  but  horizontal  winds 
measured  several  m/s.  It  is  quite  conceivable  that  a  westerly  wind  under  just  the  nght  circumstances  could 
completely  cancel  out  the  effects  of  ram  on  the  Doppler  signature,  and  eliminate  any  detecuon  of  rain  on  the 
path.  Figure  3.5  shows  the  Doppler  (hourly  median  values)  signature  for  the  month  of  March.  A  second 
plot  shows  an  estimate  of  the  impact  wind  might  have,  from  ground  measurements  of  the  wind  at  Concord. 
NH  [NOAA  Mar.  1991).  The  slow,  more  or  less  continuous,  changes  of  the  Doppler  due  to  changes  in  the 
regional  wind,  are  interspersed  with  positive  rain  signatures  showing  sudden  increase  in  the  Doppler. 
Correlauoii  ber  veen  the  graphs  is  expected  to  be  very  rough  since  the  Concord  reading  is  a  point 
measurement  at  a  different  alutude.  Nevertheless,  the  impact  of  the  wind  is  clearly  important,  especially 
during  strong  winds. 

tMore  precise  interpretation  of  the  data  may  be  achieved  by  compirison  of  the  data  with  MIT  RADAR 
measurements  made  along  pan  of  the  Prospect.  Hill  -  MLTug  path.  These  are  only  made  during  conditions 
of  rain,  however.  It  may  be  necessary  to  adjust  all  the  data,  once  a  determination  of  the  impact  of  die 
misalignment  has  been  made  -  including  the  data  from  previous  years  -  if  the  impact  is  considered 
significant. 

Another  minor  detail  that  should  be  clarified  is  that  the  PRN  code  generator  operated  at  less  than  an 
optimum  level.  Rather  than  a  30  dB  difference  (=  1023)  between  signal  and  noise,  it  has  in  reality  been 
approximately  18  dB.  (This  is  the  case  because  of  hardware  inaccuracies  in  the  state-of-the-art  PRN  code 
generators  operating  at  extremely  high  rates.)  This  reduces  the  sensitivity  of  the  system,  and  more  of  ihe 
delay  spread  measurements  during  low  signal  reception  were  discarded.  (When  the  signal  is  very  low.  the 
noise  floor  of  the  PRN  generator  would  be  ol  the  same  level  as  the  signal,  yielding  a  fiat  response  • 
indicating  an  infinite  delay  spread.) 
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Chapter  4 

Data  Collection  and  Processing  Procedure 

4 .  t  General 

After  the  initial  hardware  processing,  the  data  were  collected  on  a  personal  computer  system  on  MtTug.  It 
was  then  transferred  to  the  Thayer  "Tl"  VAX  where  it  was  permanently  stored  on  magnetic  tapes.  Most  of 
the  subsequent  processing  was  done  on  that  machine;  some  of  the  final  arrangement  and  plotnng  of  the  data 
was  then  done  on  the  Northstar  system. 

As  is  common  practice  for  troposcatter  systems,  the  data  were  analysed  in  terms  of  median  values.  It  is 
expected  that  this  will  remove  statistical  anomalities  similar  to  shot  noise,  such  as  occasional  scatter  from 
aircraft  -  see  Figure  4.1  for  a  dmesenes  of  the  mean  Doppler  shift,  showing  both  aircraft  scatter  and  a  more 
extended  period  of  rain.  Averaging  processes  are  also  expected  to  efficiently  remove  the  Rayleigh  fast  fading 
process  effects.  With  large  amounts  of  data,  the  statistics  of  the  median  value  should  approach  that  of  the 
mean  of  the  process. 
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Figure  4.1;  Timeseries  of  C  and  Ku  band  mean  Doppler  shifts.  Note  the  close  agreement 
at  the  rwo  frequencies,  except  during  rain,  where  the  Ku  band  Doppler  is  higher. 

Each  marie  numbers  the  16  measurements  made  during  an  hour. 

4.2  Data  Collection 

The  data  collection  consisted,  mainly,  of  recording  the  ui-pliase  (i)  and  quadrature  fq)  components  of  the 
power  output  of  a  120  Hz  low  pass  filter,  and  from  them  calcuLaung  the  total  power  and  Doppicr  frequenc. 
shift.  This  was  first  done  at  every  tap'  of  the  1023  long  PRN  code.  Once  this  had  been  accomplished,  a 


'window'  of  241  taps  length  was  centered  around  the  tap  with  the  maximum  power  received.  To  quicken  the 
process,  subsequent  recording  of  die  data  were  only  made  on  the  window-.  This  length  was  considered 
suf  ficient  to  aquire  an  accurate  esumale  of  die  delay  profile.  Due  to  the  relative,  unpredictable,  dnft  of  the 
codes  (caused  by  relative  drift  of  the  5  MHz  Cesium  beam  frequency  standards),  the  window  was  re<entered 
once  every  hour.  Thus  the  power  for  the  entire  1023  taps  code  was  recorded  at  the  beginning  of  the  hour, 
followed  by  sixteen  241  tap  windows  of  recorded  power  and  Doppler.  Averaging  of  the  i'  and  q' 
components  was  also  done;  60  measurements  of  each  were  averaged,  sampled  at  a  rate  of  250  Hz.  At  the 
end  of  the  hour  data  were  recorded  and  processed  for  a  spectrum;  the  processing  was  done  on  site  to  reduce 
data.  The  Doppler  and  power  measurements  were  calculated  on  site  from  the  raw  data  for  the  same  reason. 
The  data  for  the  spectra  were  collected  at  the  maximum  power"  tap.  The  data  from  the  spectra  yield 
independent  estimates  of  the  Doppler  spread  and  the  Doppler  offset.  Note,  however,  that  these  spectral 
estimates  are  relevant  to  one  tap  only,  not  to  the  enure  volume. 

There  were  few  serious  problems  with  the  data  collection  other  than  what  was  caused  by  the  relauve 
frequency  drift  A  few,  easily  recognized,  isolated  incidents  did  occur  when  faults  in  the  code  clocking  at  the 
receiver  end.  caused  the  PRN  generator  to  skip  a  random  number  of  taps.  This  had  been  a  considerable 
problem  during  the  fall,  but  was  fixed  while  the  transmitter  was  down  during  that  time. 

4.3  Data  Processing 

The  hourly  sixteen  241  taps  long  windows  formed  the  basis  for  the  stausucal  estimates  of  power.  Doppler 
and  delay  sptead.  The  total  power  is  the  result  of  summing  the  power  (subtracting  the  noise  power)  in  ihe 
window.  The  mean  Doppler  is  the  power- weighted  average  of  the  Doppler  at  each  tap.  An  esumate  for  the 
'2<j'  delay  spread  was  initially  based  on  the  variance  of  the  power  with  regard  to  the  ups.  This  was  soon 
found  to  be  insuffiently  accurate  (because  of  the  large  size  of  the  window)  and  the  distance  between  the  half¬ 
power  points  were  found  exactly.  Finding  the  peak  power  (which  was  then  halved  before  finding  the 
crossing  poults’)  involved  some  averaging.  The  esumate  for  the  delay  spread  is  therefore  a  slightly 
conservative  one.  No  estimate  for  the  delay  spread  was  made  when  the  power  was  too  low  relauve  to  the 
noise  power. 

The  noise  power  was  initially  estimated  at  the  beginning  of  every  hour;  it  is  the  lowesr  short  (5  taps) 
segment  of  powers  found  when  scanning  all  the  taps.  This  lowest  power  is,  generally,  located  far  away 
from  the  region  containing  the  signal.  This  same  criterion  was  used  subsequently  for  the  smaller  windows 
For  a  delay  profile  with  a  large  spread,  however,  no  distancing  from  the  region  containing  the  signal  is 
possible.  For  that  reason,  if  ihc  .."use  power  found  in  the  1023  long  code  sequence  was  lower  than  it.ai 
found  in  any  tap  in  the  window,  then  the  noise  power  from  the  long  sequence  was  subsuiuied  for  the  one 
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found  in  the  window.  This  reveals  an  inherent  weakness  in  the  results,  albeit  small,  since  channel 
characteristics  can  change  considerably  over  the  span  of  tens  of  minutes. 

Finding  the  total  power  and  mean  Doppler  required  finding  the  onset'  tap,  that  is,  the  tap  which  corresponds 
to  the  shortest  path  possible  across  the  mountains.  Furthermore,  no  meteorological  processes  can  be 
studied  without  an  exact  knowledge  of  the  location  of  the  onset  tap.  Finding  the  onset  is  difficult  only 
because  of  the  unpredictable  nature  of  the  relative  frequency  drift.  After  trying  various  algorithms  for 
finding  the  onset,  a  method  whereby  the  user  inputs  the  hourly  onset  peaks  was  determined  to  be  the  only 
pracucal  solution.  The  method  consists  of  the  following.  The  first  delay  profile  window  of  every  hour  in 
any  given  day  is  plotted,  as  in  Figure  3.1.  The  onset  tap  is  guessed  at  and  entered  into  the  computer.  The 
processing  program  interpolates  between  the  points  to  find  the  onset  taps  for  die  remaining  15  windows  for 
that  hour.  It  does  so  by  using  the  relevant  maximum  power  taps  and  the  user-provided  onset  taps. 


Figure  4.2:  Delay  profile  generated  by  plotting  the  points  above  the  3-dB  power  levels. 
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For  as  short  a  ume-span  as  an  hour,  die  relauve  drift  was  assumed,  for  simplicity,  to  be  linear.  The 
computer  then  plots  a  picture  for  the  whole  day,  showing  all  ‘he  points  of  die  delay  profile  above  ihc  3-dB 
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points,  as  shown  in  Figure  4.2.  In  other  words,  24  (hours)  umes  16  delay  profiles  are  plotted,  effectively 
showing  motion  and  thickness  of  significant  events.  Note  that  there  is  no  absolute  power  reference  in  this 
plot  at  any  one  ume,  points  above  the  3-dB  point  refer  to  the  maximum  found  in  that  particular  delay 
profile.  If  the  transitions  on  this  gTaph  are  not  smooth,  new  possible  onset  taps  are  guessed  at,  and  a  new 
graph  generated.  This  process  is  iterated  until  a  smooth  plot  has  been  generated.  It  is  this  process,  then, 
that  evaluates  the  validity  of  the  guesses  of  the  onset.  Once  this  stage  has  been  accomplished  can  the  onset 
taps  be  used  in  the  other  algorithms.  Exact  knowledge  of  the  onset  is  especially  important  for  the  delay 
spread  estimate,  which  otherwise  would  be  tainted  by  side-band  and  code  repetition  ’noise". 

4.4  Data  Classification 

For  ease  in  comparing  the  data  with  previous  measurements,  the  same  classifications  and  critena  were  used 
for  the  wide  band  system;  any  further  correction  of  i he  data  other  than  what  is  presented  here,  if  deemed 
necessary,  will  be  made  in  conjunction  with  corrections  of  all  prev.ous  data 

Once  the  median  values  for  each  hour  had  been  found,  the  data  were  classified  into  three  categories  according 
to  the  mode  of  propagation: 

1)  Turbulent  volume  scatter.  This  was  the  default  mode,  considered  clear  air. 

2)  Rain  scatter.  Since  the  Ku  band  is  more  senstive  to  rain,  it  was  used  to  detect  iL  If  the  median 
level  was  above  +2.5  Hz,  and  the  maximum  above  +6  Hz.  that  hour  was  designated  a  rain 
penod. 

3)  Ducting.  This  was  only  detected  when  the  C  band  was  in  use,  because  of  its  higher  sensitivity 
tc  this  phenomena.  It  was  identified  by  the  presence  of  very  strong  signals,  and  near-zero 
Doppler  shifts. 

The  data  were  then  Sectioned  according  to  hour  of  the  day,  the  quarter  day  interval  (morning:  04:00-10  00. 
afternoon:  10:00-16:00,  evening:  16:00-22:00,  night  22:00-04:00),  and  according  to  season  (Winter: 
January-Febniary,  Spring:  March-May,  Summer  June).  Note  that  EST  was  used  throughout  the 
experiment  Further,  variations  in  the  data  about  their  median  values  were  recorded  as  the  number  of  hours 
in  strong  attenuation  (>l0dB  less  than  the  median  signal  level),  attenuation  (5-10  dB  less  than  the  median  i. 
normal,  amplification  (5-l0dB  stronger  signal  than  the  median),  and  strong  amplificaaon  (  >10  d3  above 
the  median  level)  conditions.  Dunng  times  when  both  the  C  and  Ku  band  data  were  available,  a  5x5  maim 
for  these  categories  were  created  for  comparison  of  the  two  bands  dunng  different  propagation  modes. 
Cumulative  distributions  were  created  from  the  data,  for  invest  '.gallon  of  their  statistical  behaviour.  The 
next  chapter  details  these  results  and  plots  comparisons  of  the  different  quanuues  when  appropriate,  along 
with  a  discussion  and  conclusion  of  the  same. 
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Chapter  5 

Results  and  Analysis 

5.1  Introduction 

Figures  5.1  through  5.6  show  an  example  of  the  raw'  median  value  time  senes,  for  the  month  of  June. 

The  distributions  of  the  measurements  shown  is  typical  for  ail  the  months  observed.  Atypical,  however,  is 
the  extended  period  of  low  signal  received  during  the  15th  to  17t*i  of  June  -  most  clearly  seen  in  the  Ku 
band  data.  Figure  5.1-  and  in  the  extended  penod  of  ducting  during  the  19th  to  20th  of  June.  Since  this 
may  be  atypical  behaviour  for  the  channel,  and  because  of  the  limited  amount  of  data  (this  month 
constitutes  the  entire  summer  data  set),  the  variation  of  the  signal  behaviour  for  this  month  is  expected  to 
be  higher  than  normal. 

Rain  is  idenufied  on  the  path  by  its  high  Doppler  shift  in  the  Ku  band.  Figure  5.3.  Occasionally,  it  can 
also  be  discerned  by  high  Doppler  shift  in  the  C  band  This  may  occur,  as  it  did  during  the  15th  of  June, 
when  the  Ku  band  records  such  a  low  signal  that  it  loses  the  signal  on  the  1023  PRN  code.  This  occurs 
when  the  signal  is  not  detectable,  by  a  simple  search,  above  the  self-noise  of  the  PRN  generator.  When  it 
happens  -  unusual  for  extended  tune  periods  •  no  accurate  Doppler  estimate  in  the  Ku  band  can  be  made.  (In 
fact,  the  pulse-pair  algorithm  provides  such  a  sensitive  measure  of  the  Doppler  shift  that,  if  the  appropriate 
region  of  the  PRN  code  could  be  identified,  a  good  estimate  could  often  be  achieved  even  with  the  signal 
hidden  in  the  noise.)  This  represents  a  slight  change  in  the  detection  scheme  from  previous  years,  necessary 
because  of  different  noise  levels  of  the  different  systems. 

Figure  5.3  also  reveals  a  Haw  in  the  algorithms  that  determine  whether  or  not  rain  occurs  on  the  path.  The 
flaw  originates  in  the  antenna  alignment  offset  described  in  Section  3.3.  duhng  easterly  winds.  Isolated, 
large,  positive  Doppler  shifts  are  indicative  of  periods  of  rain,  as  during  the  12th  of  June.  There  are, 
unfortunately,  times  with  slow,  smooth,  changes  in  the  Doppler  shifts  -  corresponding  to  the  slow  changes 
in  wind  direction  and  speed  (see  Section  3.3)  -  causing  Doppler  shifts  above  the  rain  detection'  threshold. 

An  examp!?  can  be  seen  during  the  18th  of  June.  Even  with  the  adjustments  of  the  algorithms  that  have 
been  used  (involving  the  maximum  Doppler  shifts  for  the  hour  -  a  median  of  the  16  measurements  must  be 
above  2 .5  Hz  and  the  maximum  for  that  hour  above  6  Hz  for  the  hour  to  be  considered  a  rain  penod).  ihese 
penods,  falsely  identified  as  rain  periods,  could  in  the  worst  case  constitute  as  much  as  10-20^  of  the 
entire  rain  data  base.  Conversely,  a  strong  westerly  wind  can  counteract  and  overcome  posiuve  rain 
Doppler.  During  such  times,  rainy  conditions  would  be  considered  clear  atr  conditions.  This  probably 
occurs  less  often,  but  an  accurate  measure  is  difficult  to  obtain.  The  detection  schemes  in  this  thesis  are, 
nonetheless,  exactly  as  during  previous  years,  except  for  the  slight  alterauon  described  above.  Tlus  provides 
continuity,  and  hence  a  basis  for  comparison. 
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Previous  estimates  of  the  multipath  delay  spread  showed  much  less  sensitivity  than  the  present  estimate. 

The  range  of  the  previous  estimate  had  been  contained  within  the  100  to  200  nanosecond  lines  for  almost 
all  the  measurements.  Figure  5.5  shows  how  drastically  the  estimate  has  improved,  although  there  is  little 
change  in  the  median  of  median  values.  The  estimate  of  the  variance  about  that  median  will  improve, 
indicated  by  its  increase  from  previous  years.  This  will  represent  a  significant  increase  in  the  percentage  of 
time  that  the  available  bandwidth  is  expected  to  be  lowered  ( and  raised,  although  dial  is  of  less  concern),  as 
compared  to  previous  estimates.  This  is  because  there  are  phenomena  (ducting,  thin  layer  scattering 
dominance)  during  which  the  real  multipath  delay  spread  is  small.  However,  the  less  sensitive  two-sigma' 
delay  spread  estimate  used  previously  would  not  reveal  such  sharply  contoured  signatures.  Note  that  several 
hours  have  been  deleted  from  the  multipath  plot  This  is  done  when  the  received  signal  levels  are  too  low 
for  an  accurate  measure  of  the  spread  to  be  obtained.  This  occured,  for  instance,  during  the  15th  to  the  Pth 
of  June. 

Ducting  on  the  path  is  discerned  best  by  the  C  band  measurements.  A  low  C  band  Doppler  median 
measurement  and  an  unusually  high  signal  level  median  (also  as  compared  with  the  Ku  band  signal  level) 
qualifies  the  hour  considered  for  the  ducting  category. 

5.2  Period  of  Study 

Table  5. 1  details  the  number  of  hours  of  data  that  were  used  in  the  statistical  summaries  of  this  chapter. 

Due  to  a  break  down  of  the  transmitter  in  the  fall  of  '90,  no  data  were  recorded  until  the  beginning  of 
January  '91 .  During  the  fust  few  days  of  transmission,  a  permanent  collection  scheme  was  developed.  Data 
from  those  days  will  not  be  included  in  the  statistical  study  since  inconsistencies  would  result.  Subsequent 
collection  was  virtually  uninterrupted  through  the  middle  of  April.  At  this  time  the  local  cable  TV 
company  •  whose  antennas  are  located  at  the  receiver  site  •  claimed  that  the  project  caused  interference  with 
several  TV  channels,  wherefore  its  operation  was  temporarily  shut  down.  The  complaint  had  followed  soon 
after  the  installation  of  an  additional  amplifying  stage  in  the  Ku  band  front  end.  (This  had  been  done  in 
preparation  for  the  multi-Grequency  experiment)  The  interference  seemed  intermittent,  probably  dependent 
on  the  transmission  conditions  of  the  TV  channel.  Variou-  ouempts  to  reduce  the  interference  were  made, 
until  the  addition  of  a  narrowband  filter  finally  cUared  any  noticeable  effects  on  the  TV  stations.  For  this 
reason,  data  collection  was  intermittent  in  Late  April  and  during  May  During  June,  more  hardware 
problems  with  the  transmitter  caused  loss  of  data  during  a  few  days.  A  final  transmitter  break  down  on  June 
21st  marks  the  end  of  the  data  collection  campaign. 


Ku  band: 

Clear  Air 
Rain 

Ductina 

Winter  91 

796 

4 

Spring  91 

1164 

134 

...  10 

Summer  91 

270 

56 

33 

Total  number 
of  hours 

2230 

194 

43 

Total 

800 

1308 

359 

2467 

C  band: 

Winter  91 

Spring  91 

Summer  91 

Total  number 
of  hours 

Clear  Air 

244 

270 

514 

Rain 

29 

56 

85 

10 

33 

43 

Total 

283 

359 

642 

Table  5.1:  Summary  of  Observation  Period  (Hours) 

5.3  Comparison  of  C»band  and  Ku-oand  Signal  Levels 

Since  dve  C  band  was  used  a  reference  standard  for  the  Ku  band,  it  «s  important  to  distinguish  dissimilar 
behaviour  of  the  «:6nals.  To  this  end,  the  5x5  signal  level  matrices  in  table  5.2  are  listed,  for  the  three 
modes  cf  propagation.  In  this  case,  it  was  convenient  to  combine  the  data  from  May  and  June  -  the  only 
months  when  both  bands  were  in  operation  -  rather  than  separate  the  data  into  the  seasons.  This  is  also  an 
appropriate  step  in  order  to  increase  the  da:*  base  far  the  comparison,  especially  so  since  the  June  data 
shows  unusually  high  meteorological  activity  (see  Section  5.1).  The  categories  of  the  attenuated  and 
amplified  sta&s  are  determined  by  comparison  of  the  data  with  the  median  levels.  For  Table  5.2.  the 
medians  used  are  those  of  the  corresponding  mod  *  of  propagation.  In  the  Figure  5.7  -  5  .9  plots,  the 
median  level  shown  is  that  of  the  clear  air  observations. 


Table  5.2:  Relative  Behaviota  of  Ku  and  C  band  Signals 


S.  Alien.  Atten^ 

Normal 

Amnl. 

_  S.  AmpL 

Ku  band 

3S.  Auen. 

0 

5 

r. 

2 

0 

Alien. 

2 

30 

46 

2 

0 

Norm. 

17 

54 

199 

51 

14 

Ampl. 

0 

1 

19 

17 

24 

S.  Amp. 

0 

0 

5 

5 

!0 

Total  number  of  hours: 

514 

Ku  band  median  value: 

-118.7  dBm; 

C  band  median  value: 

•106.3  dBm 

a)  Gear  Air 


S.  Aden:  Strongly  attenuatfd  re’ative  to  the  normal  medial  Jignal  le"el  (more  than  by  10  dB), 
aitan:  menu  tied  (5  to  10  dB),  Amp:  amplified  (J  to  10  dB),  and  S.Amp:  Jtrongiy 
amplified  (more  than  10  dB). 
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■tlifTiTnHsM 

Kubarid 

S.  Alien. 

0 

0 

2 

2 

0 

Alien. 

0 

8 

9 

1 

0 

Norm, 

0 

7 

26 

2 

n 

Ampl. 

0 

2 

11 

5 

0 

S.  Amp. 

0 

0 

7 

1 

0 

Total  number  of  hours:  85 

Ku  band  median  value:  -116.2  dBm; 

b)  Rain 

C  band  median  value: 

-104.0  dBm 

BIliETuTW 

IKToMWililiS? 

Kuband 

S.  Alien. 

0 

i 

0 

0 

0 

Anen. 

2 

4 

0 

0 

0 

Norm. 

2 

5 

21 

2 

1 

Ampl. 

0 

0 

2 

0 

2 

S.  Amp. 

0 

0 

0 

0 

1 

Total  number  of  hours:  43 

Ku  band  median  value:  -ll 

3.2  dBm; 

C  band  median  value: 

-90.2  dBm 

c)  Ducting 

Freni  Table  5.2  a)  and  Figure  5.7  it  is  evident  that  the  behaviour  of  the  two  signals  is  similar  by  far  the 
greatest  portion  of  the  signal  levels  is  in  adjacent  categories  -  the  deviations  about  the  respective  median 
levels  differ  by  less  than  10  dB.  The  majority  of  the  measurements,  55%,  however,  are  in  different 
categories,  and  22%  of  those  are  in  non  adjacent  categories.  It  is  in  place,  therefore,  to  question  both  the 
appropriateness  of  the  propagation  mode  categories  considered,  and  the  methods  of  determining  the  mode  of 
propagation  for  any  given  hour.  Some  divergent  areas,  such  as  the  region  with  attenuated  Ku  band  signal 
and  amplified  C  band  signal,  may  be  statistical  artifacts.  They  could  also  be  an  indication  that  refinement 
of  the  theory  •  by  considering  more  categories  -  is  necessary.  It  is  clear,  however,  that  the  regions  with  the 
greatest  digressions  from  the  median  levels  in  Figure  5.7  -  with  amplified  Ku  band  signals  accompanied  by 
normal  C  band  signal  levels,  and  with  strongly  amplified  Ku  band  signals  accompanied  by  strongly 
amplified  C  band  signals  •  in  fact  correspond  to  the  identified  rain  and  ducting  groups,  as  evidenced  by 
Figure*  5.8  and  5.9.  No  mechanism,  however,  explains  the  very  small  group  of  data  w  here  the  Ku  band, 
unlike  the  C  band,  is  auemmed  It  seems,  therefore,  that  any  difficulties  encountered  are  primarily  die 
result  of  inexact  categorization  methods,  rather  than  the  resrJt  of  use  of  faulty  categories.  If  this  is  the  case, 
the  Figures  indicate  tha*  it  affects  only  a  small  portion  of  tie  data  With  this  in  mind,  a  conclusion  drawn 
from  these  tables  is  that  turbulent  volume  scatter  is  the  main  mode  of  propagation  at  both  Ku  and  C  bands, 
and  that  the  categories  as  well  as  categorization  methods  are  appropriate. 


Figure  5.8: 
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Figure  5.9:  Ku  vs  C  band  signal  levels.  Ducting,  May/June  1991. 

5.4  The  effect  of  Rain  on  Ku  band  Transmission 


Figure  5.8  also  shows  the  effect  of  precipitation  on  the  received  signal  levels;  for  a  larger  portion  of  the 
time  than  during  clear  air  conditions,  the  Ku  band  signal  level  is  amplified  (the  C  band  shows  less  change). 
As  expected,  (he  Ku  band  signal,  with  a  wavelength  smaller  than  that  of  the  C  band  signal,  is  much  more 
sensitive  to  the  rain  effects  (see  Eqn.l  1).  Table  5.3  clarifies  the  relationship  between  the  distribution  of 
turbulent  volume  scatter  signal  levels  to  the  distribution  of  rain  scatter  signal  levels. 


Clear  Air 

Rain 

Rain  w /  respect  t< 

hour’/  % 

hours  /  % 

Clear  Air  Median 

Median  of  all  hours 

-121.6  dBm 

-113.5  dBm 

Strong  Attenuation  (no.  hours) 

27  /  3.4% 

0/0% 

0/0% 

Attenuation 

152/19.1% 

0/0.0% 

0  /  0.0% 

Normal 

416/52.3% 

3/75% 

0  /  0.0% 

Amplification 

139/17.5% 

1/25% 

2  /  50.0% 

Strong  Amplification 

62  /  7.8% 

0/0% 

2  /  50.0% 

Total  number  of  hours 

796/100% 

4/100% 
a)  Winter 

4/  100% 

Clear  Air 

Rain 

Rain  w /  respect  i 

hours/  % 

hours/  % 

Clear  Air  Median 

Median  of  all  hours 

-120.8  dBm 

-119.5  dBm 

Strong  Attenuation  (no.  hours) 

82  /  7.0% 

6/  4.5% 

3  /  2.2% 

Attenuation 

199/17.1% 

25  /  18.7% 

21  /  15.7% 

Normal 

642  /  55.2% 

62/46.3% 

59  /  44.0% 

Amplification 

177/15.2% 

16/11.9% 

21  /  15.7% 

Strong  Amplification 

64/5.5% 

25  /  18.7% 

30  /  22.4% 

Total  number  of  hours 

1164/  100% 

134/100% 

134  /  100% 

b)  Spring 
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Median  of  all  hours 

Clear  Air 
hours  /  % 
-119.7  dBm 

Rata 

hours/  % 
-117.0  dBm 

Rain  w /  respect  ti 
Clear  Air  Median 

Strong  Attenuation  (no.  hours) 

11  /  4.1% 

2  /  3.6% 

1  /  1.7% 

Attenuation 

52/19.3% 

12/21.4% 

6  / 10.7% 

Normal 

139/51.5% 

20/35.7% 

26  /  50.% 

Amplification 

55  /  20.4% 

16  /  28.6% 

11  /  19.6% 

Strong  Amplification 

13/4.8% 

6  /  10.7% 

12/21.4% 

Total  number  of  hours 

270  /  100% 

56/  100% 
c)  Summer 

56  /  100% 

Table  5. 3:  Distributions  about  median  levels,  Ku  band,  clear  air  and  rainy  conditions. 


The  first  column  of  rain  data  lists  the  number  of  rain  hours  with  the  various  reception  levels  (as  described  in 
Section  4.4)  with  respect  to  the  median  value  of  the  rain  hours.  The  last  column  shows  the  same  hours 
with  respect  to  the  median  of  the  clear  air  hours.  As  is  seen  in  these  three  tables,  rain  has  the  effect  cf 
increasing  the  median  signal  levels,  and  reduces  the  amount  of  low-level  fields,  skewing  the  distributions 
upward.  The  rain  data  bom  the  winter  months  are  scarce,  but  the  other  seasons  contain  sufficient  data  for 
statistical  reliability.  In  all,  only  4  of  the  2467  hours,  or  0.16%,  were  dominated  by  strong  attenuation 
caused  by  rain,  relative  to  the  clear  air  median  signal  level.  Concern  for  the  effect  rain  has  on  the  path'  is 
more  appropriate  with  regard  to  high-level  interference  Hekla.  In  fact,  plot  A17  in  Annex  A,  shows  that,  for 
the  Ku  band,  high-level  fields  are  as  common  during  rain  as  during  ducting  conditions.  It  should  be  noted, 
however,  that  the  high-level  fields  correspond  to  light  rainfall,  and  that  (he  low-level  fields  are  caused,  in 
general,  by  stronger  rainfall  on  the  path.  For  design  purposes,  therefore,  it  is  insufficient  to  consider  all 
rain  as  a  single  mode  of  propagation.  A  study  must  be  made  of  the  expected  distributions  of  rainfall  on 
any  particular  path,  prior  to  system  design. 

5.5  Diurnal  Variations 


Figures  5.10  through  5.14  show  the  diurnal  variations  of  the  quantities  in  question  for  the  spring.  A 
complete  set  for  the  other  seasons  are  listed  in  Annex  A,  Figures  1  through  8.  Figure  5.10  shows  little 
variation  in  the  Ku  band  clear  air  signal  level.  During  the  mid-day  hours  there  is  amplification  of  the  signal 
by  only  2  dB.  An  increase  during  this  time  is  expected  in  clear  air  because  of  the  increased  turbulent 
activity  when  the  sun  warms  the  earth  s  surface.  The  hot,  more  turbulent,  air  rises,  reaching  the  lower  edge 
of  the  common  volume  by  about  10  am.  There  is  a  progression  of  the  diurnal  variation  that  has  been 
obscured  by  the  averaging  of  the  data  (from  March  through  May)  in  Figure  5. 10.  Each  month  shows  an 
increasing  difference  between  the  mid-day  hours  and  the  rest  of  the  day.  During  the  summer  this  difference 
becomes  quite  large,  and  is  non-existent  during  the  winter  hours  observed.  The  diurnal  variation  of  signal 
levels  during  rain  is  greater,  but  also  the  result  of  fewer  hours  of  data.  It  is  the  result  of  a  more  variable 
meteorological  process,  and  the  received  signal  level  is  more  dependent  on  the  location  and  horizontal  extent 


48 


of  the  rain,  than  is  the  clear  air  signal  level.  This  variation  during  rain  is,  therefore,  considered  statistically 
indeterminate.  The  plot  shows  yet  another  view  of  the  superior  transmission  conditions  provided  by  the 
rain.  Ducting  occured  only  during  a  few  hours  for  the  time  period  that  the  C  band  was  operational  during 
the  spring  (May).  The  signal  levels  are  clearly  higher  even  in  the  Ku  band,  not  only  in  the  C  band  (Figure 
5.11)  where  its  effect  is  the  greatest  As  can  be  seen  horn  Figure  6  in  Annex  A  (June  has  a  greater  amount 
of  C  band  data,  and  more  measurements  during  ducting),  the  greatest  coupling  into  the  ducts  occurs  during 
the  night  and  early  morning  hours.  The  effect  on  the  Ku  band  signal  levels  (Figure  2,  Annex  A)  is  as  high 
as  a  15  dB  increase  above  the  dear  air  levels.  For  the  C  band  the  increase  is  up  to  twice  that  Figure. 

Figure  5.12  and  13  show  the  Doppler  shifts  on  the  path.  The  prevailing  westerly  winds  result  in  a  clear  air 
negati  ve  Doppler  shift  signature  of  about  3  Hz.  This  is  the  result  of  path  asymmetry,  as  described  in 
Section  3.3.  Had  the  antennas  been  aligned  over  a  path  with  a  uniform  horizon  profile,  the  Doppler  shift 
would  have  been  wen  within  1Hz.  There  is  a  decrease  of  the  average  offset  throughout  the  three  seasons, 
caused  by  less  strong  and  persistent  westerly  winds.  The  Ku  band  shows  more  sensitivity  to  motion  of  the 
scanners;  the  Doppler  shift  is  roughly  proportional  to  the  signal  frequency. 

The  delay  spread  increases  during  rain,  and  is  very  much  smaller  during  ducting.  During  the  winter  (Figure 
3.  Annex  A)  this  difference  is  the  greatest,  and  during  the  summer  the  difference  is  the  smallest.  From 
experience  with  the  hourly  delay  plots  it  is  expected  that,  if  more  precise  methods  of  discerning  rain  on  the 
path  are  used,  this  difference  will  be  augmented.  The  estimates  for  the  median  delay  spreads  during  rain  in 
this  report  are,  therefore,  expected  to  be  somewhat  low. 
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Average  diurnal  variation  in  Ku  band  signal  level,  spring,  1991. 


Figure  5.10: 
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jure  5.11:  Average  diurnal  variation  in  C  band  signal  level  spring,  1991. 
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Figure  5.12:  Average  diumai  variation  in  Ku  band  Doppler  frequency  shift,  spring,  1991 
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Figure  5.12: 


Average  diurnal  variation  in  Ku  band  delay  spread,  spring,  1991 
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Table  5.4  summarizes  the  behaviour  of  the  received  signal  levels,  as  well  as  the  Doppler  shift  and  delay 
spreads,  during  June,  1991.  Similar  summary  sheets  for  the  winter  and  spring  measurements  are  listed  in 
Annex  B.  The  first  standard  deviation  value  is  the  estimator  of  the  series  of  median  values;  the  second  value 
is  the  standard  deviation  of  the  median  of  medians  given  in  the  column  above  it  The  signals  show  clear 
variation  during  the  day;  Kit  band  signals  are  strongest  in  the  afternoon  by  6  to  8  dB.  The  C  band  shows  a 
variation  of  5  to  6  dB.  There  is  a  weak  negative  Doppler  shift  due  to  the  prevailing  westerly  winds.  The 
delay  spread  shows  very  large  standard  deviation  values  compared  to  previous  years;  this  is  the  result  of  the 
improvement  in  the  algorithm  estimating  this  quantity.  The  data  for  the  Doppler  spread  values  was  not 
divided  into  these  same  groups,  but  simply  recorded  for  each  season  as  a  whole.  Previous  experience  has 
shown  that  this  measurement  is  very  insensitive  for  the  lower  range,  probably  because  of  external  noise  (60 
Hz)  affecting  the  data.  These  values  are  listed  at  the  end  of  each  table. 


Table  5.4:  Summer  1991  -  Quarter  Day  Summaries  -  Clear  Air 


Ku 

Ku 

Ku 

C 

C 

Carrier 

Doppler 

Delay 

Camer 

Doppler 

Level/ 

Shift 

Spread 

Level/ 

Shift 

Xmission  Loss 
(dBm) 

(Hz) 

(ns) 

Xmission  Loss 
(dBm) 

(Hz) 

Night 

Hours 

72 

72 

44 

72 

72 

Average 

-122.7/167.5 

-0.5 

152.5 

-106.9/151.7 

-0.4 

Std.  Dev. 

(series  /  median) 

4.4  /  0.20 

2.6  /  0.21 

70.2  /  6.5 

6.5  /  0.33 

1.5/ 0.1 

Morning 

Hours 

69 

69 

45 

69 

69 

Average 

-120.5/165.3 

•0.5 

97.5 

-102.8/147.6 

-0.4 

Std.  Dev. 

6.2  /  0.2 

2.8/0.13 

62.5/6.1 

7.6  /  0.3 

1.5/0.11 

Aftemoson 

Hour; 

65 

Oj 

35 

65 

65 

Average 

•114./158.8 

-0.1 

140 

-101.3/146.1 

0.2 

Std.  Dev. 

6.6/0.12 

2.4/0.12 

71.6/5.0 

6.3  /  0.25 

1.6/0.18 

Evening 

Hours 

64 

64 

46 

64 

64 

Average 

-120-5/165.3 

-0.3 

142.5 

-108.0/152.8 

0.3 

Std.  Dev. 

Doppler  Spread; 

6.8  /  0.25 

Hours:  264 

2.0  /  0.28 

Median:  8.6 

79.0  /  5.0 

St. Dev:  3.3 

8.7  /  0.34 

4.0  /  0.24 

Table  5.5  reinforces  the  conclusions  about  the  effect  of  rain  or  the  transmitted  Ku  band  signal.  The  signal 
levels  are  higher,  during  almost  all  of  the  day,  than  the  dear  air  measurements.  There  is  also  a  slight 
increase  in  the  C  band  signal  level  during  summer  rain.  It  is  evident  that  the  Doppler  shift  is  the  clearest 
indicator  of  rain  or  the  path.  The  rain  Doppler  spread  is  comparable  to  the  clear  air  values  during  this 
season.  There  is  insufficient  data  to  discern  any  diurnal  variations  from  this  table. 
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Table  5.5:  Summer  1991  •  Quarter  Day  Summaries  -  Rain 


Night 


Horn 

7 

7 

5 

7 

7 

Average 

-122.7/167.5 

7.3 

167.5 

-109.6/154.4 

1.9 

Sid.  Dev. 

5.3  /  1.0 

2.5  /  1.0 

89.3  /  30.2 

4.3  /  0.70 

1.2/ 

0.28 

Morning 

Hours 

8 

8 

8 

8 

8 

Average 

-114.3/159.1 

11.1 

110.0 

-109.2/154 

4.0 

Sid.  Dev. 

5.6/  1.3 

4.4  /  0.98 

84.6  /  17.7 

3.8  /  0.69 

2.6/ 

0.84 

Afternoon 

Hours 

17 

17 

13 

17 

17 

Average 

-118.7/163.5 

.8.4 

127.5 

-1023/147.1 

3.0 

Std.  Dev. 

8.8/  1.2 

4.8/  i.l 

61.3  /  10.3 

3.4  /  0.45 

2.4/ 

0.49 

Evening 

Hours 

24 

24 

23 

24 

24 

Average 

-116.5/161.3 

10.9 

140.0 

-107.1/151.9 

23 

Std.  Dev. 

Doppler  Spread: 

6.5/0.81 

Hours:  54 

4.5  /  0  54 

Median:  8.7 

68.0  /  9.7 

SLDev:  2.2 

4.7  /  0.44 

2.0/ 

0.46 

Table  5.6  shows  [he  same  quantities  during  ducung  conditions.  As  expected,  the  signal  levels  are 
considerably  higher,  peaking  during  night  hours  in  the  C  band,  where  ducting  is  of  greater  impact.  The 
delay  spreads  are  extremely  small  compared  to  that  of  either  rain  or  turbulent  scatter.  This  is  expected,  since 
the  power  scattered  toward  the  receiver  will  come  from  (he  vertically  limited  duct  Note  that  this  result  also 
applies  to  the  Doppler  spread  measurement.  It  is  much  smaller  during  ducting,  as  one  would  expect 
because  of  the  stability  of  the  ducting  region. 

Table  3&  Summer  1991  -  Queries  Day  Summaries  •  Ducting 
Night 


Houn 

11 

11 

10 

11 

11 

Average 

-113.2/158 

0.8 

30.0 

-87.3/132.1 

0.1 

Std.  Dev. 

6.5  /  0.77 

3.2/0.11 

44.0  /  22.4 

103/1.6 

0.20/ 

0.1 

Morning 

Hours 

13 

13 

13 

13 

13 

Average 

-1123/1573 

-0.2 

82.5 

-90.0/134.8 

-0.1 

Std.  Dev. 

4.8  /  0.58 

1.1/ 0.2 

46.9/11.6 

6.5/  1.2 

0.21  / 

0.45 

Afternoon 

Hours 

6 

6 

6 

6 

6 

Avenge 

-111.4/156.2 

0.7 

42.5 

-89.2/134 

-0.2 

Std.  Dev. 

3.1/0.63 

1.8/031 

17.7/  11.7 

1.7/  1.2 

0.31  / 

0.17 

Evening 

Hours 

3 

3 

3 

3 

3 

Average 

-115.9/160.7 

-1.1 

155.0 

-93.5/138.3 

-0.3 

Std.  Dev. 

Doppler  Spread 

3.4/  1.7 

Hours:  264 

0.64  /  0.22 

Median:  16 

58.5/51.0 

SiDev:  2.6 

10.9  /  2.5 

0.12/ 

0.22 
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95%  bounds.  Since  the  major  portions  of  the  curves  lie  within  these  bounds  (90%  are  expected  to),  it  is 
evident  that  the  lognormal  model  provides  an  excellent  approximation  to  the  signal  level  distributions. 
Annex  A  lists  a  more  complete  set  of  these  distribution  curves  for  different  conditions. 


Reduced  Variate  far  a  Normal  Dutnbuoon 

Figure  5.15:  Sample  cdf  for  Ku  and  C  signal  levels,  night,  spring,  1991,  clear  air. 

-70  r 


-80 


-3-2-10  1  2  3 

Reduced  Variate  for  a  Normal  Distribution 
Sample  cdf  for  Ku  and  C  signal  levels,  morning,  spring,  1991,  clear  air. 


Figure  5.16: 
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Reduced  Variate  for  a  Normal  Distribution 


Figure  5 . 15:  Sample  cdf  for  Ku  and  C  signal  levels,  afternoon,  summer,  199 1 .  dear  air. 


Reduced  Variate  for  a  Normal  Distribution 

Figure  5.!6:  Sample  cdf  for  Ku  and  C  signal  levels,  evening,  summer,  1991,  clear  air. 
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The  CDFs  are  also  applied  to  the  other  quantities;  the  Doppler  spread  and  delay  spread  cdfs  are  given  in 
Figures  5.19  and  5.20.  The  Doppler  spread  distnbutions  show  little  difference  between  rain  and  clear  <ur,  in 
the  estimate  of  their  median  values,  but  is  substantially  lower  during  ducting.  The  range  of  the  Doppler 
spread  during  rain,  however,  is  considerably  smaller  than  the  range  during  clear  air.  Note  the  full  range  of 
the  delay  spread  estimates  in  Figure  5.18,  and  the  excellent  fit  within  the  normanl  distribution  5%  and  95% 
bounds. 
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Figure  5.19: 
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Reduced  Variate  for  a  Normal  Distribution 
Sample  cdf  for  Ku  Doppler  Spread,  all  conditions  and  hours,  summer.  1991. 
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-1  0  1 
Reduced  Variate  for  a  Normal  Distribution 
Sample  cdf  for  Ku  delay  spread,  clear  air,  all  hours,  summer,  1991. 


Figure  5.20: 
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5.6  Seasonal  Variation 

The  seasonal  variation  of  three  quantities  are  of  primary  interest  with  regard  to  seasonal  change.  Again,  a 
more  complete  set  of  CDFs,  showing  seasonal  variauon,  is  provided  in  Appendix  A.  The  fust  plot.  Figure 
5.21,  shows  the  diurnal  variation  of  the  Ku  signal  levels  for  all  three  seasons.  The  summer  shows  a  great 
hump'  during  the  warmer  hours  of  the  day,  but  its  night-time  level  is  barely  higher  than  the  level  during 
the  other  seasons.  The  second  plot.  Figure  5.22,  shows  the  greater  dominance  of  westerly  winds  during  the 
colder  seasons.  With  this  in  mind,  if  the  data  is  to  be  revised  at  a  later  stage,  these  months  should  be 
investigated  first  for  accurate  categorization  of  rain  data.  (Since  the  amount  of  offset  is  a  mild  indicator  of 
the  amount  of  rain  data  that  could,  potentially,  have  been  defined,  erroneously,  as  clear  air  data.)  The  thud 
plot  depicts  die  seasonal  change  of  the  delay  spread.  As  mentioned  earlier,  the  spread  is  as  much  as  50% 
higher  during  the  -'aimer  than  winter  hours.  This  result  is  in  spite  of  the  fact  that  ihe  lowest  average 
signals  were  observed  during  the  winter  months.  (If  the  signal  >s  too  low,  there  is  a  nsk  that  the  delay 
spread  measured  will  be  unduly  large;  it  is  clear  that  this  risk  has  been  minimized  sufficiently,  by  the 
thresholds  set  in  order  for  the  recorded  data  to  be  used  in  the  statistical  presentation.)  Another  view  of  ihe 
seasonal  changes  is  presented  in  Appendix  A,  plots  A14-16.  For  all  three  cases  of  power.  Doppler  spread, 
and  delay  spread,  higher  values  were  observed  during  the  summer  season. 
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Figure  5.21; 


Diurnal  variation  of  Ku  signal  levels  for  dur:  seasons,  clear  air. 
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5. 7  Conclusion 

The  measurements  taken  during  1991  show,  in  agreement  with  results  from  the  initial  stages  of  the  data 
collection  campaign,  that  the  signal  levels  can  be  modeled  by  the  lognormal  probability  distributions,  and 
that  the  delay  and  Doppler  spreads  can  be  modeled  by  normal  probability  distributions.  The  median  values 
have  been  shown  to  be  consistent  with  the  thin  layer,  turbulent  volume  forward  scatter  model  [Crane  1991  ]. 
It  is  therefore  possible  to  predict  both  the  high  and  low-level  fields  expected  from  the  path,  by  using  the 
distributions  and  the  median  levels  predicted  by  the  theoretical  model.  Moreover,  rain  on  the  path,  unless 
large  amounts  of  heavy  rainfall  are  expected,  will  not  adversely  affect  the  wave  propagation.  Care  must  be 
taken,  however,  to  guard  against  the  variations  in  the  delay  spread  -  this  may  in  fact  be  of  more  concern 
than  signal  level  deterioration,  and  its  seasonal  variations  may  be  of  more  concern  than  the  variations  due  to 
precipitation. 

Future  work  on  the  project  would  focus  on  a  comparison  of  the  1991  data  set  with  data  from  the  previous 
few  yean,  and  a  study  on  a  monthly,  rather  than  seasonal,  basis  with  a  complete  set  of  data;  this  would 
yield  estimates  for  the  'worst  month  characteristics.  The  effects  of  the  antennae  misalignment  must  also  be 
dealt  with,  perhaps  by  comparing  the  data  with  independent  weather  RADARs.  Although  a  fairly 
substantial  portion  of  the  data  categorized  as  rain  could  be  in  error  because  of  the  misalignment,  this  is  not 
expected  to  affect  the  conclusions  made  here. 


Signal  Ixvel(dBm)  ^  Mean  Uglier  (II/) 
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Annex  A 
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Plo*  of  Avenge  diurnal  variadon  in  Ku  bend  signal  level,  winier.  1991. 
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A2  P!o<  of  Average  diumal  variation  in  Ku  band  mean  Doppler  frequency,  winter.  1991. 
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A5  Plot  of  Average  diurnal  variation  in  Ku  bind  mao  Doppler  frequency,  spring.  1991. 
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A6  Plot  of  Average  diunul  vanauon  in  C  band  signal  ievel,  spnng.  1991. 


Plot  of  Average  diurnal  variation  in  C  band  mean  Doppler  frequency,  spring,  1991 


0  6  12  18 

Hour  of  Day 

Plot  of  Average  diurnal  variation  in  Ku  band  delay  spread,  spring,  1991. 


Reduced  Variate  of  a  Normal  Distribution 


Sample  CDF  for  Ku  and  C  signal  levels,  rain,  night,  spring.  1991  (bounds  of  clear  air). 
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Reduced  Vahate  of  a  Normal  Distribution 


Sample  CDF  for  Ku  and  C  signal  levels,  rain,  morning,  spnng,  1991. 
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Annex  B 

Summary  tables  for  spring  and  winter  1991: 

Table  Bl:  Spring  1991  -  Quarter  Day  Summaries  •  Clear  Air 


Ku 

Ku 

Ku 

C 

Carrier 

Doppler 

Delay 

Carrier 

Level/ 

Shift 

Spread 

Level/ 

Xmissioo  Loss 
(dBm) 

(Hz) 

(ns) 

Xmission  Loss 
(dBm) 

Night 

Hours 

305 

305 

222 

70 

Average 

-121.4/166.2 

-1.6 

105. 

-107.7/152.5 

SuL  Dev. 

(series  /  median) 

5.9  /  .26 

3.3/. 10 

76.6  /  2.9 

6.0  /  0.3 

Morning 

Hours 

303 

303 

223 

59 

Average 

-121.7/166.5 

-1. 

105. 

-106.9/151.7 

Std.  Dev. 

6.2  /  0.25 

3.5  /O.lt 

66.8  /  3.0 

5.5  /  0.32 

Afternoon 

Hours 

263 

263 

229 

52 

Average 

-119.9/164.7 

-1.3 

97.5 

-107/151.8 

Std.  Dev. 

7.7  /  0.34 

3.3/0.14 

67.9  /  2.9 

3.7  /  .26 

Evening 

Hours 

291 

291 

250 

63 

Average 

-118.8/163.6 

•1.8 

no. 

-107.6/152.4 

SuL  Dev. 

6.6  /  0.27 

3.3/0.11 

71.7/3. 

6.0  /  0.26 

Doppler  Spread: 

Hours:  1153 

Median:  6.7 

StDev:  2.8 

Table  B2: 

Spring  1991  - 

Quarter  Day  Summaries  •  Rain 

Night 

Hours 

36 

36 

26 

3 

Average 

-122.9/167.7 

7.8 

140.0 

-100.1/144.9 

SuL  Dev. 

1U/  1.1 

3.8  /  0.38 

76.3  /  12.1 

1.4/  1.4 

Morning 

Hours 

29 

29 

28 

4 

Average 

-116.3/161.1 

9. 

145 

-108.2/153. 

SuL  Dev. 

8.1  /  1.0 

3.6  /  0.65 

79.2/  15.14 

3.2  /  0.87 

Afternoon 

Hours 

30 

30 

28 

8 

Average 

-119.9/164.7 

7.0 

122.5 

-104.4/159  2 

SuL  Dev. 

7.  /  0.97 

4.6  /  0.63 

76.2/  11.7 

4.4  /  0.69 

Evening 

Hours 

39 

39 

35 

14 

Average 

•119.8/166.6 

8.0 

157.5 

•99.5/144.3 

SuL  Dev. 

7.1/0.89 

3.9  /  0.44 

69.3  /  9.9 

2.4  /  0.66 

Doppler  Spread: 

Hours:  124 

Median:  8.1 

SLDev:  2.6 

C 

Doppler 

Shift 

(Hz) 


70 

-0.9 

2.0/0.10 

59 

-0.4 

1.6/ 0.1 

52 

0.4 

1.1/0.16 

63 

-0.5 

1.9/0.14 


3 

0.8 

1.9/0.54 

4 

1.1 

3.2  /  1.4 

8 

2.3 

3-3  /  0.86 

14 

14 

3.2  /  0.68 


Tabic  B3: 


Spring  1991  -  Quarter  Day  Summaries  -  Ducting 


Night 


Hours 

2 

2 

2 

2 

Average 

-114.5/159.3 

-2.3 

27.5 

-101.1/145.9 

Sid.  Dev. 

0.9/  1.5 

3.3/ 3.1 

72.5  /  34.9 

7.7  /  1.2 

Morning 

Hours 

4 

4 

4 

4 

Average 

-115.3/160.1 

-0.6 

45.0 

-94.8/139.6 

Std.  Dev. 

2.2  /  0.37 

0.39/0.19 

19.3  /  7.2 

4.1  /  1.2 

Afternoon 

Hours 

3 

3 

3 

3 

Average 

-115.8/160.6 

-0.8 

90.0 

-97.0/141.8 

Sid.  Dev. 

0.58  /  0.57 

0.39  /  0.25 

53.9/16.5 

0.79/2.5 

Evening 

Hours 

1 

1 

1 

1 

Average 

-123.1/167.9 

5.0 

222.5 

-99.9/144.7 

Std.  Dev. 

0.0  /  3.7 

0.0  /  5.0 

0.0/51.3 

0.0/  1.6 

Doppler  Spread; 

Hours.  5 

Median:  8.3 

St-Dev;  4.5 

Table  B4: 

Winter  1991  • 

Quarter  Day  Summaries  -  Clear  Air 

Ku 

Ku 

Ku 

Carrier 

Doppler 

Dehry 

Level/ 

Shift 

Spread 

Xmisston  Loss  (Hz) 

(dBm) 

(ns) 

Night 

Hours 

207 

207 

159 

Avenge 

*121.8/166.6 

-4.0 

87.5 

Std.  Dev. 

(senes /median) 

6.9  /  0.34 

3.3/0.13 

63.9  /  3.0 

Hours 

183 

183 

145 

Average 

-122.1/166.9 

-3.9 

70.0 

Std.  Dev. 

6.3  /  0.4 

4.2/0.18 

572  /  3.1 

Afternoon 

Hours 

183 

185 

169 

Average 

•1 19.7/164  J 

-3.3 

72.5 

Sid.  Dev. 

7.1/  0.39 

3.5/0.15 

60.1/3.1 

Evening 

Hours 

220 

220 

179 

Average 

-122.2/167 

-3.5 

80.0 

Sid.  Dev. 

62  /  0.31 

2.3/0.16 

62.6  /  3.4 

Doppler  Spread: 

Hours:  752 

Median:  7.3 

St-Dev.  2.9 

2 

-O.t 

0.35/0.6 

4 

-0.1 

0.17/0.28 

3 

0.2 

0.29/  1.5 

1 

0.4 

0.0/ 0.5 
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Table  B5: 

Winter  1991  • 

Quarter  Day  Summaries  -  Rain 

Night 

Hours 

0 

0 

0 

Morning 

Hours 

1 

1 

1 

Average 

-115.8/160.6 

375 

185 

Std.  Dev. 

Afternoon 

Hours 

3 

3 

3 

Average 

-113.5/158.3 

3.1 

107.5 

Std.  Dev. 

2.5  /  0. 

2.7  /  0. 

80.4  /  0. 

Evening 

Hours 

0 

0 

0 

Doppler  Spread: 

Hours:  4 

Median:  5.8 

SLDev:  0.7 
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MISSION 

OF 

ROME  LAB  ORA  TOR  Y 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  ail 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Comput-  r  Science  and  Technology,  Electiomagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


